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OPTIMAL SCHEBUIING IN HYDRO- THBRI4AI PO^fflR SYSTEMS 
BY THE IfflIEOD OP LOCAL VARIATIONS 

In an interconnected power system, it is possible 
to supply a given load demand in many waya» and hence 
it is natural for the operator to look for the 'best' or 
'optimum* operating strategy. Thus optimal scheduling 
in a power system implies the determination of a strategy 
which would optimize stipulated operating criterion. A 
particular choice of this criterion in power system 
operation is the cost of the fuel consumed and the 
objective is to minimize the same. Early work in this 
area was that of Kirchmayer, who solved a purely thermal 
system problem, using the coordination equations derived 
by the Lagrange multiplier technique. In this formula- 
tion, the effects of the equality and inequality 
constraints imposed by the transmission and generating 
systems were either approximated or completely neglected. 
Carpentier in 1962 formulated this purely thermal problem 
more rigorously, taking into consideration the several 
equality and inequality constraints of the system, which 



were not considered ty Zirclimayer. SuTssequently sev^eral 
solution techniques were proposed in the literature to 

solve the problem as formulated by Garpentier. However, 
the interest in this thesis is directed towards a power 
system which consists of thermal as well as hydro 
(including pumped storage) plants. The optimal schedul- 
ing of such a system is different from a purely thermal 
system in the following aspects. No fuel cost is 
associated with hydro stations and the total amount of 
wa,ter that is available for power generation over a 
specified period is limited. Hence the solution to this 
problem at any given time consists of determination of 
a plan for withdrawal of water from the hydro reservoirs 
for power generation and determination of the correspon- 
ding thermal generations so that the total cost of fuel 
is minimized. At the same time, the total power demand 
on the system is met and the operational constraints 
such as limits on the reservoir storages, the rates of 
discharges from the reservoirs, the power generations 
at various units, the voltages and phase angles in the 
network, the reactive power generations, and the line 
flows stip-ulated by stability considerations etc. are 
satisfied. The dynamics of the hydro system makes this 
a variational problam. 

Several attempts have been made in the past to 
solve the above problem by dynamic programming, incremental 



dynamic programming, and several indirect optimization 
methods like the discrete maximtim principle or the 
continuous maximum principle., , In all these attempts, 
the dif.ficiilty encountered was ‘that of enormous compu- 
tational requirements for a problem of realistic size. 
Therefore, this investigation is aimed at finding an 
algorithm, which is simple, easy to implement, and 
requires less computational effort and storage. A direct 
search method known as the ’Method of Local Variations’ 

i 

(MLV), (due to Krylov and Chernous Ko') is found to meet 
the above requirements due to its simplicity, efficiency, 
and ease of implementation. The essential features of 
this method are a discrete-time description of the system, 
a decomposition of the combined system into hydro and 
thermal subsystems, a starting nominal trajectory of 
hydro generations, and a systematic iterative method of 
perturbing the nominal trajectory such that the total 
fuel cost is monotone decreasing from iteration to 
iteration till a satisfactory convergence to an optimal 
trajectory is obtained. At every stage of the algorithm, 
all the constraints of the problem are satisfied. 

In this thesis, two mathematical models of the 
hydro- thermal system are considered, pertaining to 
(l) a short range scheduling problem, and (2) a long 
range scheduling problem. In the short range model, 
the electrical netwo^rk of the power system is represented 
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in detail. This enables consideration of various 
engineering and operating constraints on voltages, 
reDctive powers and line flows etc. for the entire 
scheduling period. , In the long range problem, however, 
only an approximate solution is sought. Instead of 
representing the electrical network in detail, the 
transmission losses in the system are approximated 
by Kirchmas^er 's loss coefficients. No other electrical 
constraints like those on voltages, reactive powers etc 
are considered, Ofcourse in both the formulations, 
constraints on hydro storages, discharge rates and 
poxfer generations are considered. 

To illustrate the applicability and advantages 
of the MLV, first a one-hydro- one-thermal problem 
similah to the one solved by Bernholtz and Graham 
(using incremental dynamic programming method) has 
been solved. Next a more complex long range problem 
consisting of three hydro and four thermal plants has 
been solved. The computational advantages of the WLY 
approach are established. 

Next a short range problem consisting of two 
hydro and two thermaJ. plants, similar to the one 
considered by Bonaert, El-Abiad and Eoivo has been 
solved. It is shown that the MLV applied to this 
problem is superior to the method used by Bonaert etal. 



il modification of tiie basic MiY approach., vfhlch consists 
in varying the step size in the algorithm is presented and 
it is shown that this can further reduce the computational 
effort significantly. 

Finally the optimisjation problem which 'includes 
pumped storage hydro plants is also considered. It is 
shown that the MLV approach permits this situation to be 
considered without any major modifications. Numerical 
results are presented to illustrate the method. 



CHAPTER 1 


IHTRODUCTIOir 

1.1 GBNERAI. INTRODUCTION TO THE PROBLEM 

With the advent of high speed digital computers, 
it has hecome possible to conduct optimal planning and 
scheduling studies of large interconnected power systems 
in a rigorous mnner using sophisticated mathematical 
models for the various components of the power system. 

In an interconnected power system with various power 
sources, more than one operating strategy is possible to 
meet the power demand in the system and hence it is ^ 
natural for the operator to seek an optimum operating 
strategy. This is determined by optimizing a pre- 
specified operating criterion with reference, to some 
control variables in the power system. An important 
criterion in a power system is the cost of the fuel 
consumed at the generating stations and it is to be 
minimized. This problem of optimal schedtiling or 
economic load dispatch has been the concern of the 
operating engineers over the last two decades. 
Eirchmayer’s [l] early work on economic dispatching in 
a purely thermal power system stimtCLated considerable 
interest in this area. His method consists of 
essentially solving a set of coordination equations 
formulated using the Lagrange multipliers technique. 
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The solution determines the most economic generation 
levels for a set of thermal generating units of unequal 
efficiency feeding into a lossy transmission system and 
constrained to satisfy the power demand. The transmission 
loss in the system is approximately considered using the 
loss coefficients [2], It was Carpentier [ 5 ] who pointed 
out in 1961, that Zirchmayer^s formulation was incomplete 
in the sense that the equality and inequality constraints 
imposed hy the transmission and generation systems were 
either approximated or completely disregarded. Further- 
more, the reactive power injection which controls the 
voltage level in the network was omitted in this formu- 
lation, Carpentier [ 3 ] in 1962 formulated this problem 
as a nonlinear programming problem, wherein the trans- 
mission network is modelled by stipulating that the 
equality constraints defined by the power flow equations 
be satisfied. Also the voltages, phase angles, and real 
and reactive injections are required to satisfy the 
specified inequality cqnstraints. The function to be 
minimized in this case is the fuel cost at the generating 
stations in the system. Several iterative techniques 
such as SUMT and gradient techniques [ 4 »5, 6] have been 
developed in the literature to solve the above nonlinear 
programming problem. 

However, a modern power system may consist of 
a large number of thermal, conventional hydro and pumped 
storage power plants connected to various load centres 



throiogh a lossy transmission network, in important 
objective in the operation of such a power system is to 
generate and transmit power to meet the system demand at 
minimum fuel cost by an optimal mix of the various types 
of plants. The study of the problem of optimum scheduling 
of power generation at various plants in a power system is 
of paramount importance in a large country like India, 
where most state electricily boards depend for their 
power requirements on both hydro and thermal power 
generations and the power systems are going through a 
phase of rapid expansion. 

Hydro-thermal optimal scheduling problem is a 
more complex one than the purely thermal system optimi- 
zation problem discussed above since the hydro reservoir 
dynamics and the constraints on the reservoir storages 
and withdrawals of water from the reservoirs have also 
to be taken into accoxint. In this problem the cost of 
hydro power generation (fuel cost) can normally be 
neglected in comparison with the cost of thermal power 
generation. However, the water that can be used over a pre- 
scribed 'interval is limited by the capacity, and inflows 
into the reservoirs, prespecified withdrawal from the 
reservoirs for meeting agricultural, navigational and 
interstate, rei^i^rements etc. The procedure for integra- 
ting the operation of hydro and thermal generation in a 
hydro-thermal system for minimum cost of generation has 



4 


"been referred to as hydro-thermal coordination [2]. The 
solution to this economic or optimal load scheduling 
problem for hydro-thermal systems consists of determina- 
tion of a plan for withdrawal of water from the reservoirs 
for power generation, and the corresponding thermal 
generations so as to minimize the total cost of thermal 
generations over a specified period, while meeting the 
total power demand on the power network. The dynamics 
of the hydro system makes this a variational problem. 
Itirther, the resulting schedule has to satisfy the 
operational constraints such as limits on the reservoir 
storages, discharge-rate constraints, limits on voltages 
and phase angles, limits on active and reactive generations 
and limits on the power transfers along the individual 
lines etc. These constraints are primarily imposed on the 
variables in order that the equipment ratings and the 
desired operating conditions of the generators and the 
electrical network are not violated. 

Due to the exhibition of cyclic nature, by the 
inflows into the reservoirs and the load demand on the 
power system, the problem of optimal scheduling in 
hydro -thermal power systems can be broadly classified 
into long and short range problems. The long range 
problem is normally defined over an year and the short 
range problem over a day or a week. A search into the 
literature reveals that the basic difference between 
these two problems is the manner in which the 



characteristics of the power system network are represen- 
ted in the mathematical model used in the analysis. In a 
long range problem, the transmission losses in the system 
have been approximately considered as a qxaadratic function 
of the power genera.tions in the system, using the loss 
coefficients due to Zirchmayer [2]. The operating con- 
straints on voltages, reactive power generations, the line 
flows etc. in the power network are not considered. But in 
a short range problem the transmission losses in the system 
are considered by the accurate solution of the nonlinear 
power flow equations of the power network. In this 
formulation, the operational constraints on the voltages, 
reactive powers, the line flows etc. are duly incorporated. 

The proper integration of hydro and thermal 
generation in a power system of realistic size is quite 
complex and presents considerable difficulty in the 
determination of the optimal hydro and thermal schedules. 

Many algorithms have been proposed in the past for this 
problem. In this thesis another algorithm is proposed 
which the author believes requires less computational 
effort and storage than the existing ones. Iny such 
claim can of course be established only after a comparative 
study is made of the existing and the proposed algorithms 
bringing out their relative advantages and disadvantages. 

Some effort is made towards this direction in the present 
work, however it is difficult to do due to lack of sufficient 
information regarding data and numerical results in the 
literature. 
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1.2 STATE OE THE ART OE THE PROBLEM 

Severs,! attempts have been made in the past to 
solve the optimal scheduling problem for hydro- thermal 
power systems using direct as well as indirect search 
techniques such as Dynamic Programming, Pontryagin’e 
maximum principle etc. A detailed survey of the'more 
important solution methods used for this problem solution 
is given in [?]. A brief overview of the state of the 
art of the above problem is given below. 

Dandeno, Chandler and Eirchmayer [8] assumed that 
the hydro plants operate at essentially constant head and 
formulated a short range problem with loss coefficients 
as follows. Let E(P (t)) be the cost • of -■ generation 
in unit time for a generation of P|j,^ units at the o-th 
theimial generating station. Let W^(t) be the discharge 
rate at the i-th hydro station for a generation of Pg^ 
units and let t)e the total volume of water available 
for power generation. Let Eg(t) be the load demand at 
any time t and let transmission losses 

expressed as a function of thermal and hydro generations 
using loss coefficients. Then it is desired to 

tf QJ 

minimize / ( I E(P«.(t)) > dt (1.1) 

0 j=i 

with the restriction that 

¥. (t) dt = E. 

0 ^ 


f j P 


( 1 . 2 ) 



and 


|l ^ i=L 


where a = niimber of steam plants 
p-a= number of hydro plants 

t^ is the total interval of optimization*. By the applica- 
tion of variational methods, the following coordination 
equations are derived 




aPij(t) 




(1*4 ) 




( 1 . 5 ) 


where 




ap-J-t) , . . 

ID plant j 


incremental production cost of steam 


— * incremental transmission loss of steam 


Tj 


plant j 


■ g = incremental transmission loss of hydro 

plant i 


d¥..(t) 

ap^(t) 


incremental water rate at hydro plant i 


X = incremental cost of received power 

= water conversion coefficient which 
converts incremental water rate into 

plant cost 
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l=a+l 




( 1 . 6 ) 




loss coefficients. 


Bquations (1.2) to (1.5) when solved give the solution, 
for Pm.(t) and P„-(t) over the entire interval.’ The 

X 3 lix 

values which effectively convert incremental water 
rates into incremental plant costs, determine the 
amotmts of water used at each hydro plant. Therefore 
must he so chosen that the desired amounts of water 
are utilized, i.e. eqn. (1.2) is satisfied. ^ is so 
chosen that the load demand is always satisfied in 
eqn, (1.5). It is observed that s are constants and 
X is a function of time. 


When the number of hydro plants to be coordinated 
increases, a choice of Y conversion factors, which will 
result in desired water storage becomes very difficult. 
Similarly the choice of X will also be difficult as 
computational instability was observed by Bandeno [9] 
while iterating on X for the loads at certain hours. 
Bandeno [ 9 ] further observed that by direct application 
of the coordination equations, solutions are obtained 
which sometiiaes dictate generations outside the plant 
capacities, because these constraints are not included 
in the problem formulation. 



BeTnholtz and G-raliam* [lO] to [I3] coiisiiiieDred a 
short range protilem with, loss coefficients sad attempted 
to solve it tlirough Incremental Dynamic Pro^raoiiing (IDP), 
procedure. This method consists essentially 0* applying 
the Dynamic Programming algorithm locally, in. a neighhour- 
hood of a nominal trajectory. Starting from an initial 
feasible trial schedtile, the problem of deteriiiaing 
optimum schedule for a system consisting of aiQr number 
of constant head hydro stations and any number* of steam 
stations, was reduced to a sequence of prolDleiUiS involving, 
first, the determination of hydro-schedules "whiicch 
maximize the weighted output of the hydro system and second 
the determination of the corresponding mSmimuim cost 
thermal schedules. With increased nujmber of fajr-dro 
stations in the system, this method becomes mmeh involved 
as the computational requirements become enoraieasly high. 

Sokkappa [14] formulated the long raiijge problem 
(see Sec. 2.1) as a nonlinear programming prOBtlem. Trans- 
missioa losses were considered through the usse of loss 
coefficients. Por each subinterval of optiniLzstion , 
the constraint, which is most likely to be viioHated was 
picked up and a slack variable associated wltfci. this 
constraint. The gradient of the cost fum-ction. was 
evaluated and the steepest descent method, us-ei to obtain 
the solution of the problem, starting from a Icnown initial 
schedule. Such a procedure beccmes impossible in view of 
large dimensional requirements for problems of laulti- 
reaervoir hydro-thermal problems,. 
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Dynamic Programming was used ly Petersen [15] to 
solve an annual optimization problem. Hano etal [l6] 
and Dahlin and Shen [l7] have described computational 
approaches to the economic operation of a multi- 
reserVoir hydro-thermal system employing the continuous 
maximum principle of Pontryagin. Uarita [l8] and Oh [l9] 
have employed the discrete maximum principle to solve the 
same problem. In all these methods, the computational 
requirements for a realistic size problem become 
prohibitively high. 

Bonaert etal [20] have employed the method of 
solution given by Bernholtz and Graham [10| in the 
framework of a more sophisticated model. By a suitable 
decomposition of the total system into a hydro subsystem 
and a thermal subsystem, their method consists of 
starting with a feasible nominal hydro schedule, the 
thermal subsystem is solved for optimal thermal 
generations, by obtaining optimal power flow solutions 
[6] one for each subinterval. The incremental costs of 
hydro powers known as dual variables are also obtained 
from the optimal power flow solutions. Using these 
dual variables, a new neighbourhood hydro schedule is 
obtained in such a way that the overall cost function 
is minimized in each iteration. The iterative process 
is continued till the overall minimum cost of operation 
is obtained* As evidenced by the example given in [20], 
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it appears that this method also sxif f ers from the large 
computational requirements for practical size systems. 

G-opala Rao [21] has formulated the short range 
problem as an additively separable nonlinear programming 
problem and used lasdon's decomposition [.22] technique, 
to split the problem of larger dimension into subproblems 
of smaller dimensions to be solved iteratively. The 
main difficulty in this method is the initial choice of 
the dual variable vector, which plays a dominant role 
in the computational time required. 

1.3 OUTLINE OF THE THESIS 

The central theme of this work is the application 
of the Method of Local Variations [23h [24] (MLV), to 
obtain solutions to the complex problem of optimal 
scheduling in hydro-thermal power systems. A brief 
outline of the material discussed in the following 
chapters is given below. 

The second chapter introduces the two basic 
approaches which arise due to the repetitive nature 
of the system data, like the load demand on the system, 
the inflows into the reservoirs etc. This leads to the 
classification of the problem into short range and long 
range problems. Justification of the consideration of 
the transmission losses in both the cases is attempted. 
Tiio mathematical models have been formulated defining 
the problem of optimal schedtiling in hydro-thermal power 
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systems, oae for tcheasLort range problem and the other 
for the long ran gee prooboLem. The main difference in 
these t-wo noaeUs lie tJ:he manner in which the electrical 
network is repn:asse!i~1:e»sd- A brief discussion is also 
included to "br ing o'^totthe salient features for and 
against both t lie cciosd&ola described In this chapter. 

Ctiapt er* 3 Ifbr iefUy descrihes the computational 
techniciues u_sed im tbie solution procedures adopted in 
later chaptero. Ihae L’ethod of Local. Yariations in the 
context of i-ts :apopl.iocivtion to the work reported is 
described in. cle tanil. Next the Augmented Penalty 
Function apEproa-chn «of : Best enes [253 ss applied to 
mathenaticaH jr-ogrsi'iiiijBng problems by Miele etal [26] is 
discussed bn’LeQgy. (Opf;imal Power Slow solution procedure 
of Dommel ami !Iiinia_ejrv f[6] is then briefly summarized. The 
material in. fcliris cstep~ter is aimed at providing a good 
introduction to tlae 1) asic features of the solution 
techniques eapBo oyed. , 

In Ghsp ter | the long range model is considered. 
Solution to “tlaiiB aoo*£lel is obtained using the Method of 
Local Varist iena, i numerical example is given to 
illustrate tie i^Ouxitf. on procedure for a one-hydro- one- 
thenaal pr oblem. . Hi this case a nominal hydro schedule 
is asstmed to s'ta-xtftwith and the corresponding thermal 
generations aaoe csttftaiined by solving a quadratic equation 
of power balaaic ce cooiudition f or each subinterval • 



The nominal hydro sohed^io is perturbed 1 b its ueigh- 
hourhood usiBg the MhV iteratiueXy to ohtaiu the opti^. 
eoheduie. A oOBpariscB hetveeu the MLV pros edure and 
the method due to Bemholtz aud draham [lO] is sxven 

the above problem. The method has then beer 

1+4 multi thermal problem, 

extended to solve a malti hydro, mui 

Im owing the nominal hydro generations a 
In this case, Knoviag 

-1 r^+ riPter mined from the nominal trajectory, 
each hydro plant determinea 

it is necessary to solve a sequence of mathematical 
programming problems one ior each sublntervul. Both 
equality and inequality constraints have been tahen 
iato account in solving the mathematical programmrng 
problems using the Augmented Penalty Bunotlon anproao 
adopted by Hlele etal [ 26 ] and the modifications given 
in sec. 5.3. A numerical example of a power system 
consisting of three hydro and four thermal plants xs 
30 lved and the results are presented. A recommended 
procedure for the choice of the step size vect 

, . rT-i+h the MIV is made based, on the 
■while -worhing with xne 

computational experience gained during the course of 
this work. 

In Chapter 5, a short range problem is considered, 
in this case a detailed mathematioal model is employed, 
wherein the transmission losses are accounted for y 
solving the A.C. power flow equations during each soh 

X 1 Tn each subinterval, knowing the^^^ ^ 

ling subinterval. In eacn suoix 

j 4 - ciT the'hvdro plants, en optii®! 

hydro generations at all the nyar p 
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power flow solution [6] is obtained for solving the 
optimal thermal generations. While obtaining the 
optimal power flow solution, the nominal hydro generations 
are considered as injections into the power network at 
the buses at which the hydro plants are situated. Thus 
starting with a nominal trajectory of hydro schedules, 
a sequence of optimal power flow solutions are to be 
computed one for each subinterval. The nominal trajec- 
tory is then improved iteratively by perturbing it in its 
neighbourhood by using the MIV. This ensures a monotone 
decrease of the cost of the thermal generations in the 
system, from iteration to iteration. A numerical 
example illustrates the procediare. A brief comparison 
of this method with that of Bonaert etal [20] is also 
given. 

Chapter 6 shows the generality of the MLY, by 
applying it to solve a power system optimal scheduling 
problem consisting of conventional as well as pumped 
storage hydro plants. Numerical example illustrates 
the method. 

The seventh chapter gives an overall summary of 
the work reported in this thesis and a few suggestions or 
possibilities for further research in this area are given. 



CHAPTER 2 
PROBIEM DEFINITION 

2.1 INTRODUCTION 

The problem • of optimm schedioling in hydro- 
thermal power systems ean be broadly classified into 
long and short range problems. The long range problem 
refers normally to an annual problem* The inflows 
into hydro reservoirs exhibit an annual cyclicity. 
Furthermore, there may be a seasonal variation in 
power demand on the system, and this too exhibits an 
annual cyclicity. The optimization interval of' on© year 
duration is thus a natural choice for long range 
optimal generation scheduling studies. The solution 
to the scheduling problem in this case consists of 
determination of the amounts of water to be drawn from 
the reservoirs for hydro-generation in each scheduling 
subinterval and the corresponding thermal generations 
to meet the load demand over each subinterval, utilizing 
the entire amount of water available for power genera- 
tion during the total in tervsil. It is necessary in this 
problem to predict beforehand the load demand, inflows 
into and evaporations from the hydro-reservoirs for the 
entire optimization year. In this case the load demand 
on the system is probabilistic in nature. For the 
purpose of scheduling studies normally forecasts are 
prepared from the past load data of the system using 
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statistical prediction methods. Reservoir dynamics 
should he considered in the formulation of the problem 
[27] since there may be significant variation in 
reservoir head in the optimization period. The normal 
scheduling subintervals for this problem are a day» a 
week or even a month and average values for hydro and 
thermal generations in these subintervals are determined 
in the long range problem. Due to very nature of the 
problem rigorous power network model representation is 
not necessary in this case. It is sufficient to 
consider the transmission losses in the system using 
loss coefficients developed by Eirchmayer [ 2 ]. 

The short range problem usually has an optimisa- 
tion interval of a day or a week. This period is 
normally subdivided into hourly sub intervals for 
scheduling purposes. The solution to the long range 
problem will stipulate the amounts of water to be 
utilized over each day or week (the scheduling sub- 
interval of the long range problem). Then the 
solution to the short range problem consists of an 
optimal plan for utilization of this water for power 
generation and the corresponding optimal theimal 
generations determined, considering the load demand 
on the system and the constraints imposed on its 
operation. In this case, it is the normal practice to 
assume that the load demand and the inflow data for the 
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reservoirs are known deterministically. Approximate 
determination of transmission los|!es "by utilizing the 
loss coefficients, as in done in 4 long range scheduling 
problem, may not yield schedules sufficiently close to 
the optimal ones. Furthermore the loss coefficients [2] 
are determined for a given system network configuration 
as a function of the generated powers in the network. 

Also there is no way of checking the operating limits 
on the voltage magnitudes in the network and the 
reactive pox^er generations at their sources. However 
for the practical implementation of the optimum schedule, 
it is necessary to incorporate the constraints on the 
electrical network such as the upper and lower limits 
on the bus voltages, the maximum and minimum allowable 
reactive power generations by the reactive power sources, 
and the limits imposed on the phase angle difference 
between ajiy two adjacent buses for stability considera- 
tions. Thus the transmission losses are therefore 
usually included exactly by solving the nonlinear power 
network equations. The head variations in reservoirs 
during the optimization interval of a short range problem 
are usually insignificant and can be neglected. There 
may however be situations where this variation has to 
be taken into account by considering the reservoir^^^^ 
dynamics. The above two models of the optimal hydro- 
thermal scheduling problem are f ormulated .in Secs . 2. 2 
and 2.5* 



2.2 MATHEMATICAL MODEL POR LOKG RAKGB SCHBDULIITG 


A power system consisting of *h’ hydro plants 
and ’r* thermal plants is considered. A discrete-time 
representation of the continuous time optimal schediil.ing 
problem is presented here. The total Interval of optimi- 
zation is divided into N equal suhintervalsr ' subinterval 
being assumed to be each of unit length for simplicity. 

The hydro stations are assigned to be operating with 
reservoirs which are independent of each other. It is 
also assumed that the load demand Rjj(k), reservoir 
inflows Lj^(k) and evaporation losses ej^(k), (k=l,...,N), 
over the optimization interval are known accurately enough 
in advance. The initial and final reservoir storages 
represented by the vectors X(l) and X(N+l) respectively 
are specified. These values represent the constraint on 
the amount of water to be consumed over the optimization 
interval.The state equations representing the hydro system 
dynamics are 


XjL(k+l) s 3Cj^(k:) + L^(k) - u^(k) - e^(k), (i*l,i ,h) , 

(k = 1,...,N) (2.1) 


The optimization problem is to determine the rates of 
discharges u^(k), (k = 1, . . . ,11) , (i = l,...,h) and the 
thermal generations P^^(k), (k = 1, . . . fH) , (i = 1» • * • »t) 
such that the performance index, representing the total 
thermal power generation cost. 


f f p^(p..(k)) 

k=l 


( 2 . 2 ) 



19 


is minimized s'ubject to the eq'uality constraint 

P:jj(1c) = ■*“ * ^Iioss^^^ ' 

(k=l,...,N) (2>5) 

and the following inequality constraints 


u, * ^‘ix.(l£)iu# p (i^lf««#3^h)^ (k 

X min ^ X X max 


^ (k) X. rn ^ f (l = 1> ♦ # • >lx) | (k 

1 min i 1 max 


min- max * * l,...,r). 


(2.4) 

2,. .. ,N) 

(2.5) 


(k = 


( 2 . 6 ) 


The power halance equation (2*5) states that the total power 
generated should be equal to the sum of total load demand 
and transmission losses in any subinterval k. Equations (2.4) 
to (2.6) specify the upper and lower bounds on the rates of 
water discharge, water storage, and thermal generations, in 
each subiiiterval. In eqn. (2.2), the function yields 
the fuel cost for thermal generation of the j-th station 
and is considered to be a quadratic of the form 




(2.7) 


where a . and b . are constants determined from the fuel 
3 d 

cost-output characteristics of the j-th thermal station. 

In eqn, (2.5)» Ej^Qgg(k), the transmission loss during the k-th 
subintertral, is expressed, following Eirchmayer [2j, in 
terms of hydro and thermal power generations as 



P. (k) ^ Y"" P,(k) B_P.(k), (k = (2.8) 

loss ^^2, j=l ^^3 3 


where 


P, (k) = P„, (k), (i = l,...,h) , 
X ill 

A 


P^(k) = (i = h+1,. ,h+r), (a = (2.9) 

The loss coefficients. . are assumed to he available for 

-.31 j 

the given power system. The hydro generation at the i-th 
station during the k~th sub interval can be related to the 
water head, the storage of the reservoir at the beginning 
and end of the interval, and the discharge rate by a 
relation [l6] of the form 


H 


Pgi(k) = --|^[1 + -~(x^(k) + x^(k+l))]u^(k), 
(i “ l,...,h), (k ~ 1,«..,B^} 


c. 


( 2 - 10 ) 


where c^^ and G are known constants. 

Thus the problem is to determine the drawdown at 
each hydro station and the generation at each thermal 
station over the optimization interval such that the water 
available for hydro generation is fully utilized, the 
equality and inequality constraints specified in eqns.(2.3) 
and {2,4) to (2.6) are satisfied and the fuel cost in 
thermal generation is minimized. 

2.3 MATHSMiTICAI- MODEL POR SHOET RANGE SCHEDULING 

As in Section 2.2, a discrete version of the 
continuous time optimal scheduling problem is employed in 
the formulation. An electric power system of *n’ buses 



21 


with ’h' hydro stations and *r' thermal stations is consider- 
ed, .33ie total interval of optimization is divided into N 
equal suh intervals, each of which is assumed to he of unit 
length for simplicity. The load in each suhinterval is 
assumed to remain constant at all load buses. It is also 
assumed here, without loss of generality, that the power at 
a specific bus is generated either by a thermal or by a 
hydro plant, but not by both. The hydro stations are 
assumed to be operating with reservoirs which are independent 
of each other. The initial and final storages of the 
reservoirs are specified. These values, which specify the 
amount of water available for power generation during the in- 
terval of optimization, are generally obtained from the 
solution of the long range problem. 

¥e will assume here for the sake of generality that 
the reseirvoir head variations are not hjegligible. Thus the 
hydro system equations remain the same as in Sec. 2.2 and 
are repeated here for the sake of continuity. State 
equations describing the hydro system dynamics therefore, 
are 

s:^(k+l) = x^(k) + L^(k) - Uj_(k), 

(i=l,...,h), (k*l,...,N) (2.11) 

Here the evaporation losses term is omitted, since in 
short range problem the evaporation effect is usually 
negligible. The hydro generation during the k-th 
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sutinterval at the i-th station can be related to the water 
head, discharge rate and the snbinterval terminal storages by 
a relation [l6] of the form 

c. 

= -^[1 + + x^(k+l))]u^(k) , 

(i = l,...,h), (k = l,...,lT) (2.12) 

Bomds on the volximes and the rates of discharge of 
water from the reservoirs are given by 


^i min“ inax ’ - 1» • • • »li) , (k - 1,,..,U) 


X. . < X. (k) < X. 

1 min 1 — 1 max 


(2.13) 

f (i “ l,**.,h), (k — 2,»«»,1T) 

(2.14) 


The formulation of the thermal system and the electrical 
network follows closely the Carpentier-Sirouxs' formxilation 
as given by Peschon [28], During any subinterval k, the 
solution to the thermal system and the electrical network 
can be viewed as a static optimization problem under the 
assumption that the generations and the loads during this 
sub interval are constant. 


The cost of power • generation in the system is 
considered to depend only on the active powers Pj^(k), 
generated by thermal plants^ Thus the performance index J, 
representing the total power generation cost 
¥ r 

j = I I P (P„ (k)) (2.15) 

k=l i=l ^ 

is to be minimized subject to the following equality and 
ineqmlity constraints. 
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Equality Constrainta: 

Equality constraints arise due to the requirement of 
real and reactive power "balance at each "bus and for each 
subinterval* Thus for a typical bus i, we have the following 
equations. 

The real power balance is given by 


I. (k) - P. (k) + a, (k) = 0 , (1 = 1 n), (k=l,...,K) 

^ ^ ^ ( 2 . 16 ) 


and the reactive power balance is given by 
Zi(k) - Qi(k) + D^(k) = 0 , (i = l,...,n), 

(k=l,...,lT) 

where 

n 

i.(k) = I v.(k)Tr (k)T aos(e + 6 (k) - 6 (k)) 

- 

+ 7^(k) T^^Cos (k * 1,...,E) 


( 2 . 17 ) 


(2,18) 


n 


E^(k) 


I ? {k)y (k)T Sin(6 + 6 (k) 
1 3 IJ i 

oVi 


6 .(k)) 


+ v?(k)Y..Sin 0.. , (k = 1,...,I) (2.19) 

1 11 11 

It may be observed that to get the expressions for and 

K., the model for the transmission line is taken to be an 
1 

equivalent % network. 

Inequality Constraints s 

The related variables Tv(k), Q^(k)> and T^(k) 
are subjected to upper and lower limits » which are 
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primarily determined "by the equipment ratings and operating 
conditions of the generators and the electrical network. 


P?(k) + C^(k) - S, 


2 < 
i max — 


P. . - P^ (k) 
1 min i^ 


%. mm" f 


Qlft) - Qi 

min- i 

man- 


0 

0 

0 

0 

0 

0 


( 2 . 20 ) 

(2.21) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


The restrictions on the power transfers along the lines 
for ensuring system stability can be represented by 

\d . (k) - 6 .(k)| < I T. . I (2.26) 

1 0 10 niax' 

Thus the problem consists of minimizing (2.15) subject to 
eqtiality constraints (2.16) and (2.17) and inequality 
constraints (2.20) to (2.26). Also the constraints imposed 
on each reservoir storage and discharge rates (2.12) and 
(2.13) have to be met. 

2.4 DISCUSSION 

The first model defined in Sec. 2. 2 relates to a 
long range problem and ija this model the transmission losses 
are approximately considered using the loss coefficients [2] • 
Further, there is no reflection of the characteristics of 
the electrical network in this model and hence there is no 



provision to adhere to the operating constraints on the 
bus voltages, reactive powers, and the line flows. The 
model is simple and the computational requirements in 
solving the problem will be quite less. But the solution 
on implementation may result in inadmissible voltage 
solutions at the buses in the transmission system and 
stability requirements stipulated in terms of the maximum 
phase difference across a particular line may not be met. 
However, this model is very useful to solve a long range 
problem which provides a good starting data base for the 
short range problem, which uses a rigorous and sophisticated 
model as defined in Sec. 2. 3* 

The model defined in Sec. 2.3 is applicable for a 
short range problem. Since in this formulation, the limits 
on the bus voltages and phase angles, reactive powers, and 
line flows are considered, the solution obtained will 
always be admissible* However, for this rigorous and 
sophisticated model, the solution methods are more complex 
and require large computational effort. 

2.5 SUMTiAHT 

In this chapter, the problem of optimal scheduling 
in hydro-thermal power systems has been mathematically 
defined following the two basic approaches of considering 
the characteristics of the electrical network. In ^ 

first model, which is normally used for the solution of 
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a long range prololem, "the transmission losses are approxi- 
mately considered, but the characteristics and constraints 
on the transmission network are totally ignored. Thus the 
solution obtained may not be always feasible. The second 
model is formulated for a short range problem. The 
transmission losses are exactly included accounting for both 
active and reactive components. The various operating 
limits on the bus voltages, reactive powers, and the line 
flow requirements are considered in this case. The solution 
thus obtained is always feasible and implementable. A 
brief comparison of the above two models is given from 
their applicability and limitations considerations* 



CHAPTER 5 


SOLUTION TECHNIQUES 

3.1 INTRODUCTION 

To solve tHe optimal laydro- thermal scheduling problem, 
the Method of Local Variations [23]> [243»^ the Augmented 
Penalty Eunction Method [26] and the Optimal Power Plow 
Solution [6], have been employed after suitable modifica- 
tions. The algorithms as used in the solution are briefly 
discussed below. 

3.2 METHOD OP LOCAL VARIATIONS (MLV) 

The Method of Local Variations (MLV) was developed 
by Chernous* Ko [23] and Krylov and Chernous’ Eo [24] for 
the numerical solution of optimal control problems. The 
method was also used by Singaraj [29] in solving optimal 
structural design problems. The major advantage of this 
method is its simplicity and its efficacy to deal with the 
various bo-undary conditions and inequality constraints on 
both the state and control variables. The essential features 
of this method are a discrete time description of the system, 
a starting nominal trajectory, and a systematic iterative 
method of perturbing the nominal trajectory, such that the 
cost index (objective function )is monotone decreasing in^^^^^^ 
successive iterations till a satisfactory convergence to 
the optimal trajectory is obtained. 



statement of -the Problem: 

Consider the discrete time version of a continuous 
time control process described by the system of difference 
equations 

2(k+l) - g(X(k),u(k),k), (k = 1,,..,N) (5*l) 

where X(k), and u(k) are state and control vectors of 
dimensions h and m respectively. N is the number of 
discrete time eubinteivals into which the tbtal interval of 
optimization is divided. It is required to find a phase 
trajectory z^(k), (i = 1,. ..,h), (k = 1,...,N+1) and its 
control u^(k), (i = l,,..,m), (k = 1,...,^), such that 
(5-1) is satisfied and the functional 

N 

J(X, u) = I PT,(X(k), u(k), k) (5-2) 

k=l ^ 

is minimized where Pj^(X(k), u(k), k) is the cost incxirred 
in the k-th subinterve.1* subject to the constraints on the 
state and control variables given by 

u(k) e U(k)» (k=l,,..,N) (5-3) 

X(k) e G(k), (k * 1,...,N+1) (3-4) 

where U’(k) and G-(k) are generally variable closed regions 
of h and m-dimensional spaces respectively. The initial 
and final conditions for the vector ^(k) are taken care of 
by properly defining the regions G-(l) and G(N+l). !H3te 
solution of the above problem Using the Method of Local 
Variations will be clear from the algorithm described below. 



29 


Algorithm Descriptionj 

An initial nominal trajectory (k - lr**»»^+l)» 

is considered to start with. This initial trajectory need 
not he a feasible trajectory. Peasihility here means that 
the nominal trajectory does not violate the state constraints 
and the corresponding control variables also do not violate 
the constraints imposed on them. A siirvey of the initial 
trajectory at the instants k s= 1,...,N+1 , reveals at 
what instants the constraints (3»3) and/or (3*4) are not 
satisfied. If for some k, (l£ k< N+l)i either (3»3) ot 
( 3 . 4 ) are not satisfied, the initial nominal trajectory is 
not feasible and hence is not suitable. In this case, the 
initial nominal trajectory is modified (sub controlled) by 
changing the unsatisfactory X^(k). This change consists 
of adding one after another to the unsatisfactory components 
z^Ck), of the vector ^(k), a quantity where the 

numbers j, r take on the values j = l,...,h, r = ^l,...,+s. 
Ax. > 0 are the given magnitudes of the steps in the 

u 

coordinates and s is a given integer. If this method of 
variation permits revision of ^(k), which satisfies the 
constraints (3*3) and (3-4), the survey of the successive 
points is carried out. If the initial nominal trajectory 
cannot be corrected by this variational (subcontrol) 
operation, at any instant, the initial trajectory is 
discarded and a different initial approximation is 
selected, which may yield a feasible nominal trajectory 
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to start with. Here the values of the cost functional 
J^(^, u^) corresponding to ^(k), (k = l,...,ljr+l) as well as 
t^DS sub interval costs are computed and stored. 

After obtaining an initial nominal trajectory, which 
is feasible, in the second part of the algorithm, the current 
nominal state trajectory ^(k), (k = 1,...,H+1) in the i-th 
iteration is varied as described below to obtain the nominal 
trajectory for the (i+l)-th iteration. The variation is 
carried out successively for each of the components of 
X^(k); initially for the first one, then .for the second 
etc. For each component x^(k), its own step of variation 
Ax^(k) > 0 is chosen. Let the initial and final conditions 
on the states i.e., X(l) and X(1T4-1) be assumed to be known 
or specified. Then it is to be noted that the initial and 
final values of the nominal state trajectories always 
satisfy the specified conditions. How the following local 
variations operation is performed on the nominal state 
trajectory at time instants k = • let it be 

assumed that variation of the state values at the k-th 
instant is under consideration, and also let x^(k) be under 
consideration.' The upper index gives the iteration number 
and the lower index gives the number of the component of 
the state vector. Completion of the variation of all the 
components of the state vector X(k), (k = 2,...,N), 
signifies the end of an iteration. At the beginning of 
the variation of the j-th component of the vector X(k), 



all the preceding components have already been varied; 


i.e., the state vector at the k-th instant is 

x^(k),...,xj(k))^. Here 

I *1 

T(k) represent the nevr values assigned to the 

J •“ J- 

first (j-1) components of the state vector in the i-th 
iteration hy the application of Itihe MLV-* Variation of the 
j-th component of the state vector at instants I- 1,.. .,k~l» 


have already been performed so that the current state 
trajectory is (r = l,...,j-l), ( S- = 2,.*.»N); 

(^= 2,...,k-l), (t =k,...,N); 

V J ^ 

(r = j+1, . . - ,h) , ( = 2, . . . ,N) . Now the value of the 

component state at the k-th instant is varied to 
iL 

2:--(k) + Ax^- (k) ,keeping all the other state values unchanged, 
lests are performed to see whether (1) this state value is 
admissible, (2) the control inputs required to translate 
the state from X(k-l) to Z(k+l) via Z(k) are admissible 
(where the states X(k-l), Z(k) and Z{k+l) are defined in 
the above trajectory), (3) the changed performance index value 
for- the subintervals (k-l) and k taken together is less 


than the corresponding performance index value for the same two 
sub-intervals for the trajectory, before the j-th component 
at the k-th instant is varied. If the answer is in the 
affinnative for the above three tests, then the j-th 


component of the state vector at the k-th instant is 

changed to x^'^^(k) = x^(k) +Ax^(k). If the answer is in 
»3 J V 

the negative, the same three tests are performed for 
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.x^(k) ~ Ax.(k). If this proves successftil the j-th 
J J 

component of the state at the k~th instant is changed to 

x^'^^(k) = z^(k) “Ax.(k). Otherwise the old value of the 

component of the state x^(k) is retained at this instant 

and the local variation operation is performed on x.(k+l) 

J 

at the instant (k+l), and so on. When the variation of 
the j“th component at all instants k = is 

completed, the variation process is shifted to (j+l)-th 
component of the state and so on. In this way, the local 
variation operation is applied successively for all compo- 
nents of the state vector at all instants of time to 
complete one iteration and obtain the nominal trajectory 
for the (i+l)-th iteration. These iterations ^ which 
ensure a decrease in the values of the cost function in 
each iteration, are performed till a satisfactory convergence 
to an optim\M has been obtained. The resulting X (k), 

(k = l,w..,]I+l) and u*(k), (k = 1 ,...,N) are the desired 
optimal state and control vectors respectively. 

3.3 AUGMENTED PENALTY FUNCTION METHOD 

Considerable work has been done over the past 
several years on the solution of minimization problem 
with equality constraints, which may be stated as 

minimize f = f(x) ( 3 * 6 ) 

subject to the constraints 


^ (x ) = 0 


(3.7) 



Here, f is a scalar function of x, x an n-dimensional 
vector, is a q-dimensional vector witli q < n. It is 
assumed that the first and second partial derivatives of 
f and <j) with respect to x exist and are continuous. 

Ariy Xq which satisfies constraints (5*7), is 
said to he a feasible solution and the set of all the 
feasible solutions is called a feasible set, 

I = { 2 ^ ^(x) = 0} (3»8) 

X is said to be an optimal solution to the problem if the 
following is true 

f(s*) — f(x) for every x e T (3«9) 

Ihe usual method of solution for the above nature of 
problem has been the method of Lagrange Multipliers. 
Hestenes [25] modified it by augmenting the function by 
a penalty function in addition to the usual terms obtained 
by the product of Lagrange Miiltipliers and the equality 
constraints (Augmented Penalty Function Method). It has 
been explored by Miele etal [ 26 ] in connection with the 
ordinary first order gradient algorithm. The method usually 
exhibits faster convergence than the ordinary Lagrange 
Multiplier method [26].- A brief description of the steps 
in the algorithm is given below. 

From the theory of maxima and minima, the problaa 
defined by eqns . (3» 6 ) and (3* 8) can be recast as that of 
minimizing the augmented .function 

^'(x, 1) = f (x) + ( 3.103 



y-T 


Here is a q—vector Lagrange mul'biplier and liie 
superscript T denotes the transpose of a vector or a 
matrix. If 

- (3.11) 

denotes the gradient of the augmented function, the 

optimum solution (z*,^ ) must satisfy the simultaneous 
equations 

i (x*) = 0 (3.12) 

K*) = 0 (3.13) 

Eqn.(3.13) is known as the optimum condition. How in 
the Augmented Penalty Function Method, a new function 
W(x, 1 , it) is defined as follows 

(3.14) 

(3.15) 

Here the scalar H > 0 is the penalty constant. It was 
sxiggested hy Miele etal [26], that faster convergence is 
possible if the Lagrange multiplier and the penalty 
constant are chosen at each iteration so as to obtain 
certain desirable properties, for instance (a) descent 
property on the augmented penalty fvinction, (b) descent 
property on the augmented function, (c) descent property 
on the constraint error, and either (d ) constraint 
satisfaction on the average or (e) individual constraint 
satisfaction.Properties (d) and (e) are employed to 
first order only* 


^(x, A » H) = 1'(2»A ) + ^(x) + E l(x) 

A {p 

where P(x) = (x) <{>(x) 
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The Lagrange mxxltiplier ^ is determined by 
minimizing the error in the optimum condition defined in 
eqn.(3»l6) with respect to _Xfor given Owing to the 
fact that the error in the optimum condition is 

Q(£,i ) = [f^(2) + $|(£) + «|(s)x ] (3.16) 

the Lagrange multiplier vector is derermined by the 
relation 

Qx(2, X ) = 0 (3.17) 

which yields 

4>y(x) X+ <l>^(x) ^(x) * 0 (5-18) 

This linear vector equation is equivalent to q linear 
scalar relations in each of which the only unknown is the 
Lagrange m-ultiplier and can easily be determined. Similarly 
the penalty constant R is chosen in such a way that on the 
average, the constraints are satisfied to first order [ 30 ] 
and is given by 

R « 2P(2)y^^(x) (3.19) 

Thus the Lagrange multiplier X^ and the penalty 
constant R are determined in each iteration as explained 
above. The basic algorithm is summarized as follows. 

The Basic Algorithm: 

Let ^ denote the nominal point, x, the vaxied 
point and A x the displacement leading from the nominal 
point to the varied point. Let A. denote the Lagrange 
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multiplier, E, the penalty constant, and a the gradient 
step sise. The gradient algorithm Is represented hy 


S' (s, i) = 

E^(S) = 2^(5) i 

W (x, X , R) = Sx^-’ i. ) + ^ 

^ -a 


( 3 . 20 ) 


X = X + 


Por a given nominal poini x» Lagrange multiplier 
and penalty constant R, equations (3*20) constitute a 
complete iteration leading to the varied point x for a 
specified gradient step size a. 


In the above algorithm it is assumed that x can 
assume any value. But if in addition to the equality 
constraints given hy eqn.(3.7), there are inequality 
constraints on x of the foim 


Then the Ineguality constraints egp-O.Zl) can he handled 
easily hy assuring that the adjustment algorithm in (5,20) 
does not send any component of x beyond its permissible 
limits. If the oorreotionAxj_ from (5-20) would cause 
to exceed its limits, set to the corresponding 

limit. 


( 3 . 22 ) 



^i max ' 

if 


X^> 

^i max 

^i = 

^i min ' 

if 


X. < 
X 

^i min 


U* * 

otherwise 




Even when a component of x has reached its limit, its 
component in the gradient vector must still be computed 
in the following cycles because it might eventually 
back off from the limit. 

In the present work, while working with the long 
range problem discussed in Sec. 2. 2, it is required to 
solve a sequence of mathematical programming problems one 
for each subinterval to solve for the optimal thermal 
generations, knowing the hydro generations from the 
nominal trajectory. In each sub interval k, the cost 
of the thermal generation (2.2) is to be minimized 
subject to the equality constraint (2.3) > (power balance 
equation) and the inequality constraints (2.6). The 
above method of Augmented Penalty Eunction is employed to 
solve the mathematical programming problem encoTintered in 
each subinterval (Ref. Chapter 4). 

3.4 OPTIMAL POWER PLOW SOLUTIOR 

In Sec. 2.5» the model for short range optimization 
problem in hydro-thermal power systems was described. It 
can be seen that the problem in each subinterval once again 
consists of minimizing a cost function ( 2 . 15 ) subject to 
the equality constraints (2.16) and (2.17) (power flow 
equations) and the inequality constraints (2.20) to (2.26), 



It can be seen that for a power system of n buses, usually 
a set of 2(n-l) equality constraints will result- We 
could employ the method of Augmented Penalty Function 
discussed in the previous section also for this problem- 
However, the Optimal Power Plow Solution due to Pommel 
and Tinney [6] which is based on simple Lagrange Multiplier 
Method would be more suitable here and will become clear from 
the description given below. 

Equations (2.16) and (2.17) generally will give 
rise to a set of 2(n-l) equations corresponding to an 
n-bus system (the equations for the slack bus would not 
be needed). Any solution which satisfies the above 
equations constitutes a power flow solution. Dommel and 
linney [6] solve the power flow problem by Newton’s 
method and then the adjustment of the control parameters 
so as to minimize the objective function is done by the 
gradient method. The optimization problem developed 
depends on the assumption that initially the inequality 
limits are not incorporated* The problem of optimal 
power flow can thus be stated as 

minimize f(x, u) (5.23) 

1^1 

where f(x, u) is the objective function and z and u 
are dependent and control variable vectors respectively.^^^^^^^^^^^^^^^^^^^^^^ 
Any bus in a power system network is characterized by 
four quantities, namely the magnitude of the voltage at 
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the bus, its corresponding phase angle, the real and 
readtiYe powers entering the bus of which two of them 
are specified and the other two are to be found. This 
will give rise to the dependent and control ■ variables ♦ 
The minimization in eqi!i.(3.25) is subject to equality 
constraints 

:s) = 0 (3«24) 

Using the Lagrangtan formulation# the above problem is 
equivalent to minimizing the unconstrained lagrangian 

I(x,u) = f(x,u) + g(x,u) (3*25) 


\ is the Lagrange multiplier vector. The following 
necessary conditions are to be satisfied for a minimum. 


3x 3£ 


+ 


3u “ 

t = 


i rn 


3g 
3 u 


T 


= 0 



( 3 * 26 ) 


(3.27) 


(3.28) 


3 g 

it is to be noted that — is the Jacobian matrix in 

3 X 

a power flow solution by Newton’s method [ 3 I] and the 
satisfaction of eqn. (3.28) is Simply an ordinary power 
flow. Hence if a power flow program by Newton’s method 
is available, can be readily obtained , Thus by 
solving the three nonlinear equations (3.26), (3»27) 
and (3.28) iteratively, the minimum can be obtained. 
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using eqn. (3*27) for updating the control variables. The 
iterative steps in the solution algorithm using the 
gradient method are as follows 

1. Assume a set of contro parameters u 

2. Find a feasible power flow solution by Newton’s 
method 

3* Solve from (3*26) for \ , 


4* Insert 2 !l fi'om (3*29) into (3.27) and compute the 
gradient vector 


Vf =: X (5.30) 

- 3u 8u3u- 

5. If is sufficiently small, the minimum has 

been reached. 

6. Otherwise find a new set of control parameters from 


^ ^oia _ ^ 


(3-51) 


and return to step 2. 

c in ( 3 . 31 ) is the step length in the gradient 
direction. 

A careful observation of eq.n.(3.27) reveals that 
it holds only if no inequality limit is imposed. Thus, 
at the optimum where some of the variables are on their 
limits, ^ 0, The handling of inequality limits 

for control variables offers no difficulty as they can 
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be set and held to their limits upon violation. When the 
dependent variables have violated their limits, the penalty 
method is used to force the variables back to their limits. 
The reason for using this approach are that these limits 
are not usually hard limits and that penalty terms can 
easily be incorporated. 

This method offers a .very neat and efficient 
algorithm for solution of the particular type of problem 
discussed with a .caref-ul choice of the step length c and 
the value of the weighting factor for penalty functions. 

5.5 SUMH4RT 

A. discxission of the solution techniques namely, 
the Method of Local Variations with suitable modifications, 
Augmented Penalty Function Method and the Optimal Power 
Flow Solution is provided in this chapter. 



CHAPTER 4 

SOLUTION TO THE lORG RANGE SCHEDUIIIG PROBLEM 

4.1 INTRODUCTION 

In this chapter, the solution to the long range 
prohlem of optimal scheduling- in hydro-thermal power 
systems is obtained using the Method of Local Variations 
and the Augmented Penalty Function Method, already 
discussed in Chapter 5. The mathematical model used 
here has been described in Sec,. 2^2^. Where the trans- 
mission losses have been approximately considered 
using the loss coefficients and the limits on the 
operating voltage magnitudes, reactive powers and line 
flows etc, are not considered. 

As already pointed out in Chapter 1, several 
attempts have been made in the past to solve the above 
problem using the method of Lagrange Multipliers, 
Incremental Dynamic Programming ( IDP ), Dynamic 
Programming, and the Pontryagin's maximum principle 
(continuous as well as discrete). The method of 
Lagrange multipliers used by Dandeno etal [8] does 
not consider the limits on the generations and hence, 
the solution obtained may not be practicable. Further 
the choice of the Lagrange multipliers and the 
Y-conversion factors (see Sec, 1.2) becomes very 
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difficult as the number of hydro stations increases. 

The formulation of the problem by Dynamic Programming 
suffers from the excessive computer storage require- 
ments even for smaller systems. The formulation by 
the Pontryagin*s maximum principle (both discrete as 
well as continuous versions) is essentially a control 
problem and becomes quite mwieldy for large systems. 
Pur t her in all the indirect methods as this, the 
constraints on the state variables present consider- 
able difficulty in handling them. Out of the various 
attempts made in the past,. the IDP procedure adopted 
by Bernholtz and Graham [10] has been found to be most 
promising and hence a direct comparison of this method 
with the proposed method in this chapter is attempted 
in Sec. 4« 6, 

The solution to the above problem as given here 
essentially consists of a discrete time description of 
the system, a decomposition of the combined system into 
a hydro and thermal subsystems, a starting nominal 
trajectory of hydro reservoir storages and a systematic 
iterative method of perturbing the nominal trajectory 
such that the cost of thermal generations is monotone 
decreasing from iteration to iteration till a satisfac- 
tory convergence to an optimal trajectory is obtained. 
At every stage of the algorithm, all the constraints of 
the problem are satisfied. In each iteration, starting 



with a nominal trajectory of reservoir storages^ 
knowing the predicted load demand and inflow data, 
the discharge rates of the reservoirs during the 
suh interval as well as the corresponding hydro 
generations are computed. Having fixed the hydro 
generations over the entire interval of optimization, 
it remains to solve a sequence of mathematical program- 
ming problems (one for each subinterval) to obtain the 
corresponding thermal generations for each subinterval 
during the minimum cost of thermal generation. Thus 
this problem becomes a two level optimization problem 
and the details are disci;issed in subsequent sections i 

For purposes of illustration (i) a one hydro 
and one thermal problem and (ii) a three hydro and 
four thermal problan have been consideredi A detailed 
comparison of the MLV with IDP used by Bernholtz and 
Graham is provided and a suitable recommendation of the 
choice of the step size in the MIV algorithms is also 
made based on the computational experience. 

4.2 PORMUIATIOH OP A OHE-HIDRO-OHE-THERMAI PROBLEM 

Consider an optimal scheduling problem of a one 
hydro and one thermal system. Let the total interval 
of optimization be divided into N equal subintervals 
(however in the method considered, the subinteryals 
need not be of equal length). Let the hydro station 
state equation be 

x(k+l) = x(k) + L(k) - u(k) -■e(k), (k=l,.,.,H) (4.1) 



x(l) and xdJ+l) are specified from the requirement of 
the storages in the reservoirs at the beginning and 
end of the scheduling interval respectively, fiie 
performance index 

k=l ■ 

is to be minimized subject to the equality constraint 
P^(k) - Pjj(k) - Pqi(I^) + ^loss^^^ " (4.3) 


and the inequality constraints 

x^i^ < x(k) < x^ax ’ 

^min- ^niax ' 


(4.4) 

(4.5) 

(4.6) 


(4.7) 


"•min max 

p„ . < Prn(l^) - » (k=l,...,l^) 

•^Tmir ^ 

where 

PCP^Cl!:)) = aria's:) + WjCls:) 

where a and t are oonstanto to te detemined from 
the fuel coat-output charaoteristios of the thermal 

station. 

P (k) = f I P,(k) B.f ^(k), (k=l,...,N) 

- iil jii i. ^ 


(4. 8) 


and ^ 

phC’^) = V ^ 
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4,3 DESCRIPTION OP THE SOLUTION ALGORITHM 

The algorithm starts with the choice of an 
initial state trajectory of reservoir storages, 
x^(k), (k=l, , . . ,N+l) . The choice of initial trajectory 
can he made keeping in view the physical considerations 
(note that the trajectory need not he a feasible one 
to start with). Prom the assmed nominal trajectory, 
the associated discharge rates and the corresponding 
hydro generations in each suhinterval are calculated 
using eqns.(4.l) and (4.9) respectively. The constraints 
on the storages and the corresponding discharge rates 
are verified. Having determined Pjj(k), (k=l,...,N), 
eqn,(4,3) becomes in this case, a quadratic equation 
in P^(k) for each k. It is solved for each k to obtain 
the corresponding thermal generation P^(k). A choice 
has to be so made in this case from the two solutions 
of the quadratic equation that the upper and lower 
limits on the thermal generation are satisfied and the 
one which gives lesser cost for the subinterval under 
consideration. The sub interval costs as well as the 
total cost are computed and stored. If at any stage, 
the limits on either the reservoir storage, discharge 
rate, or the thermal generation are violated, the 
algorithm corrects the initial assumed trajectory by 
changing the unsatisfactory state variables suitably 
as described in Sec, 3.2 until a feasible trajectory 


is obtained. 
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Now the initial nominal trajectory (feasible) 
is varied iteratively nsing the following local 
variations procedure to obtain the optimal trajectory. 
Consider the variation of the nominal trajectory 
x^(k), (k=l, . . . ,N+1) in the j-th iteration to obtain 
the nominal trajectory (k=l,^, . . ,K+1) for the 

(j+l)~th iteration. Since the initial and final 
reservoir storages are specified, the local variations 
operation is to be performed at time instants i=2,...,l!r. 
let it be assumed that variation of the state value at 
the k-th instant is under consideration. Variation of 
the state values at time instants i = 2,,,.,k-l have 
already been performed, so that the current state 
trajectory at instants 2 ,35r is x^(i), (i=2, . . . ,k-l) , 
x^(i), (i=k, ...yN). Here x^(i), (i=2, . ,k-l) represents 
new values assigned to the nominal state trajectory in 
the j-th iteration, by the application of the MIV at 
time instants i = 2,...,k-l, Now the value of the state 
at instant k is varied to x^(k) + Ax, keeping all other 
state values unchanged. Tests are performed to see 
whether (l) this state value is admissible (2) the 
discharge rates required to take the state x^(k-l) 
to x^(k+l) via x^(k) + Ax are admissible and (3) the 
changed performance index value for the two sub intervals 
(k-l) and k taken together is less than the oorresponding 
performance index value for the same two subintervals 
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for the trajectory x^(k-l) ->• x^(k) x^(k+l) If the 

ansvjer is in the affirmative for the above three tests,,, 
then the state value at the k-th instant is changed to 
x^(k) = x^(k) + A X. If the answer is in the negative, 
the same three tests are performed for x^*(k) -A x* If 
this proves successful, the state, value at the k-th 
instant is changed to x^(k) = x^(k) - Ax, Otherwise 
the old va.lu.e is retained at the k-th instant and- the 
local variations operation is performed at instant 
(k+l). In this way the Local Variations operation is 
applied successively upto the time instant N, to 
complete the iteration and obtain the nominal trajec- 
tory for the (j+l)-th iteration. These iterations 
which ensure a decrease in the values of the functional 
in each iteration are performed till a satisfactory 
convergence to an optimum has been obtained. 

4,4 VARIATION OP THE STEP SIZE IN THE METHOD OP 
LOCAL VARIATIONS ALGORITHM 

In the MLV algorithm discussed in Sec, 4.3» 
it is assumed that the step size Ax (i.e,, the step 
size used in the variation of the hydro trajectory) 
is held constant throughout the iterative process. 

One would tend to feel that the choice of Ax has to- 
be based on the engineering judgement taking into 
consideration the desired accuracy in the final^^ ^ 
schedule. However, during the course of this work 
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it was Disserved that, starting with a very large step 
size ( Ax can be as large as the maxim'um allowable 
storage of the reservoir itself) and gradually reducing 
it towards the optimim resulted in a fast convergence 
in all the cases studied. The procedure using a 
variable step size is as follows. Convergence is 
obtained initially using a broad neighbourhood around 
the nominal trajectory, that is one corresponding to 
a large step size Ax and then considering the 
trajectory thixs obtained as the nnmmni.1 trajectory,, 
the variation process is continued employing MLV to 
obtain another better approximation to the optimxmi 
trajectory with a smaller neighbourhood (say with a 
step size equal to half the previous value) and so on, 
until a satisfactory convergence to the optimum is 
obtained. This approach yielded the optimum trajec- 
tory with significantly less computational time than 
that corresponding to a small constant step size used 
throughout the iterative process, 

4.5 HIMERICAL EXAMPLE 

In this section, a numerical example which is 
similar to the one solved by Bernholtz and Graham 
is solved using the MLV algorithm des crib ed above for 
a one-hydro-ohe-thermal problem. The results will be 
presented for both constant step size and variable 
step size procedures as indicated above. The problem 
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considered infact is a short range scheduling problem 
of one day subdivided into 24 hourly sub intervals. 
Bernholtz and Grahams’ [lO] formulation corresponds to 
the approximate model discussed in this chapter. The 
motivation for this is that a comparison of the MIV and 
the IDP procedures can be made (Sec. 4.6). The data as 
taken from Ref. [lO] is given below. 


Initial storage = 

Pinal storage = 

Maximum a.llowable storage = 

Minimum allowable storage = 

Maximum allowable discharge 
discharge = 

Minimum allowable 
discharge = 

Optimization interval 
considered = 


711,000 ft ^/sec.-hour 
0.0, ft ^/sec.-hour 

711.000 ft ^/sec.-hour 
0,0 ft^/sec,-hour 

56.000 ft^/sec. 

0,0 ft^/sec, 

24 hours (one day) divided 
into 24 hourly sub intervals 


Fuel cost function used for the thermal station is 

f(P^) = 373.704 + 9 . 6 O 644 P 5 . + O.OOI 991 IP 5 millions of 
BTu/hour. Cost of fuel is assumed to be 35 cents per 
million BTu, 

P(P^) = 0,35 PtPg,) dollars 

loss coefficients : B^^, ~ 0.00015, B22 = 0*00005, 

®12 0 . 00001 . 
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Maximum 

thermal 

generation = 

400 

MW 

Minimum 

thermal 

generation = 

0.0 

MW 

Maximum 

hydro 

generation = 

450 

MW 

Minimum 

hydro 

generation = 

0.0 

MW 

Basic head 


=0 

134. 

0 feet 




c ■= 

2.0xl0“'^..._ 




G = 

0,845x10"^ 


Inflo-ws into the reservoir and evaporation from the 
reservoir are not considered in this problem, fable 4,1 
gives the initial nominal trajectory of the hydro reservoir 
storages. The optimal schedules obtained for this problem 
using the constant step size as well as variable step size 
procedures indicated in Secs. 4.3 and 4.4 respectively are 
illustrated here. Table 4.2 shows the optimal schedule of 
discharge rates, hydro and thermal generations with a 
constant step size of Ax = 500.0 ft ^/sec.-hour. Table 4.3 
gives the optimal schedule of discharge rates, hydro and 


thermal generations using the variable step size procedure, 
starting step size being A z = 711*000,0 ft^/sec,— hour , 

Table 4.4 provides the optimal trajectories corresponding 
to constant step size and variable step size procedures. 


Computations are perfonaed on IBM 7044 Computer at X.I.T, 
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TABIE 4.1; IlTITIAl NOMINAL TRAJECTORY OF HYDRO RESERVOIR 
STORAGES. 


S.No. 

Instant k (Beginning 
of the honr) 

Initial storage in 
ft ^/sec.>?*hour 

1 

1 

711,000 

2 

2 

688,000 

3 

3 

653,000 

4 

4 

620,000 

5 

5 

590,000 

6 

6 

568,000 

7 

7 

536,000 

8 

8 

502,000 

9 

9 

470,000 

10 

10 

436,000 

11 

11 

403,000 

12 

12 

376,000 

13 

13 

343,000 

14 

14 

310,000 

15 

15 

278,000 

16 

16 

245,000 

17 

17 

212,000 

18 

18 

180,000 

19 

19 

148,000 

20 

20 

116,000 

21 

21 

90,000 

22 

22 

68,000 

23 

23 

47,000 

24 

24 

23,000 

25 

25 

(end of 24th hour) 

0 
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TABLE 4-2: OPTIMAL SCHEDULE FOR A OIB-HIDRO-OHB-THEEMAL 

PROBLEM (for a constant step sizeAx=500,0 ft^/ 
sec. -hour) 


Hour 

Load 

Demand 

station generation 
in MW 

Optimal discharge 
rate in ft^/sec. 


in M¥ 

Hydro 

Thermal 


1 

460 

341.93 

137.49 

26,500 

2 

430 

321,12 

125.95 - 

25,000 

5 

420 

313.32 

122.93 

24,500 

4 

410 

311.97 

113-99 

24,500 

5 

400 

304.28 

110.90 

24,000 

6 

410 

309.27 

116.47 

24,500 

7 

470 

345.53 

144.42 

27,500 

8 

550 

393.68 

182.67 

31,500 

9 

660 

440.98 

253,65 

35,500 

10 

670 

438.12 

266.56 

35,500 

11 

680 

435.27 

279.49 

35,500 

12 

700 

432.42 

302.84 

35,500 

13 

580 

351.12 

252.33 

29,000 

14 

600 

361.23 

263.73 

30,000 

15 

610 

365.16 

270.48 

30,500 

16 

610 

363.05 

272.41 

30,500 

17 

700 

419.91 

314.11 

35,500 

18 

740 

417.06 

358.44 

35,500 

19 

700 

390.95 

340.44 

33,500 

20 

690 

382.63 

337.62 

33,000 

21 

640 

345,70 

319.54 

30,000 

22 

600 

315.10 

306.42 

27,500 

23 

550 

279.29 

288.18 

24,500 

24 

500 

243.97 

269.92 

21,500 

Initial cost 

in dollars 

,= 23/727.0 


Final 

cost 

in dollars 

= 20,399.6 
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TiBIE 4.5: OPTmL SCHEDULE EOR A OHE-HIDRO-^HE-THBEMAI. 

PROBIBM (for a variable step size starting 
with Az = TlljfOOO ft^/sec.-ho-ur) 


Hour 

Load 

Demand 

station generation 
in M¥ 

Optimal discharge 
rate in ft^/sec. 


in M¥ 

Hydro 

Thermal 


1 

460 

355.17 

125.43 

27,528 

2 

430 

328.90 

118,82 

25,612 

5 

420 

315.38 

121.04 

24,668 

4 

410 

302,55 

122.86 

23,751 

5 

400 

287.71 

126,23 

22,694 

6 

410 

298.81 

126.14 

23,668 

7 

470 

331.35 

157.42 

26,362 

8 

550 

382.76 

192.55 

30,611 

9 

660 

439.85 

254.65 

35,389 

10 

670 

441.80 

263.27 

35,778 

11 

680 

433.45 

281.12 

35,332 

12 

700 

436.05 

299.58 

35,777 

15 

580 

357.69 

246.30 

29,528 

14 

600 

368.73 

256.86 

30,612 

15 

610 

370.25 

265.81 

30,917 

16 

610 

368.09 

267.79 

30,917 

17 

700 

419.65 

314.35 

35,471 

18 

740 

416.80 

358.68 

35,472 

19 

700 

389.72 

341.56 

33,389 

20 

690 

382.05 

338.15 

32,943 

21 

640 

349.59 

315.94 

30,332 

22 

600 

320.22 

301.65 

27,943 

25 

550 

281.51 

286.08 

24,694 

24 

500 

245.25 

268,72 

21,612 

Initial cost 

in dollars ;= 

25,727. 

0 

Pinal 

cost 

in dollars =: 

20,595. 

5 • 
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TABLE 4.4; OPTIMAL TRAJBCTORIES OP RESERVOIR STORAGES 


S.Ho. 

Instant k 
(Beginning 
of the hour) 

Optimal storage 
in ft^/sec.- 
honr nsing a 
constant step 
size A3&=500,0 
ft^/sec.-hour 

Optimal storage 
in ft^/sec.— 
hour using a 
variable step 
size a2c= 711,000 
ft^/sec.-hoTxr 

1 

1 

711,000 

711,000 

2 

2 

684,500 

683,472 

5 

5 

659,500 

657,860 

4 

4 

635,000 

635,192 

5 

5 

610,500 

609,441 

6 

6 

586,500 

586,747 

7 

7 

562,000 

563,079 

8 

8 

554,500 

536,717 

9 

9 

503,000 

506,106 

10 

10 

467,500 

470,717 

11 

11 

432,000 

434,939 

12 

12 

396,500 

399,607 

15 

15 

361,000 

363,830 

14 

14 

332,000 

334,302 

15 

15 

302,000 

303,690 

16 

16 

271,500 

272,773 

17 

17 

241,000 

241,856 

18 

18 

205,500 

206,385 

19 

19 

170,000 

170,913 

20 

20 

136,500 

137,524 

21 

21 

105,500 

104,581 

22 

22 

73,500 

74,249 

25 

23 

46,000 


24 

24 

21|500 

21,612 

25 

25 

(end of 24^ 
hour) 

0 

0 ■: 



In addition to the results given for the constant 
and variable step sizes in Tables 4*2 and 4.3» two more 
trials with one constant and one variable step sizes are 
obtained and the computational time required for these 
are summarized below in Table 4.5* 


TABLE 4.5: BBBEGT OB STEP SIZE ON THE COB-IPUTATIONAL TIME 
(one-hydro-one-thermal problem) 


Step size ax 
in ft^/sec,- 
hour 

Cost in dollars 

Computational 
time in secs. 
(Execution) 

Initial 

Final 

711,000+ 

25,727.0 

20,395,5 

15.75 

10,000"^ 

23,727.0 

20,397.0 

17.77 

500 

23,727.0 

20,399*6 

20.45 

250 

23,727.0 

20,380,6 

45.76 


Note; + in the above table Indicates a variable step size, 

4.6 COMPARISON OF THE INCREMENTAL DYNAMIC PROGRAMMINC PRO- 
CEDURE WITH THE METHOD OF LOCAL VARIATIONS FOR 
A ONE-HYDRO-ONE-THERMAL PROBLEM 

Bellman* s dynamic programming algorithm proved to 
be inefficient for practical size problems, mainly because 
of the enormous computational requirements especially in 
terms of fast access memory in computers* Many modifica- 
tions of the basic dynamic programming algorithm requiring 
significantly less storage have been developed in the past. 
Incremental Dynamic Programming (IDP) was one such method.. 
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developed and successfully applied to the scheduling prohlem 
under consideration by Bernholtz and Graham [10], It consists 
essentially of repetitively applying the dynamic programming 
algorithm locally, in a neighbourhood of a nominal trajectory, 
The -word incremental refers to the fact that in each iteration, 
the search for an improved hydro schedule is confined to a 
neii^bourhood of the nominal hydro schedule obtained in the 
preceding iteration. 

The Method of Local Variations is yet another direct 
method employing search in a nei^bourhood of a nominal 
trajectory as in Bernholtz and Grahams’ method and having 
less computational requirements than the latter. Thus it 
can be seen that both of these methods involve direct search 
in a neighbourhood of a nominal trajectory. It is therefore 
natural to attempt a comparison of these methods from the 
aspect of computational requirements. In IDP, the best 
trajectory in a defined neighbourhood of a nominal trajectory 
is obtained by applying Bellman’s principle of optimality. 

The resulting trajectory is made the nominal trajectory for 
the next iteration. This process is repeated to obtain 
convergence to a (locally) optimal trajectory. The MLV 
tries to obtain in each iteration only a better trajectory 
than the existing nominal trajectory in the neighbourhood 
under consideration. 



a :- — — - — — i.> . tiJiie 

a b c d 

i X 

f h 

PIOURE 4.1: NOMINAL TRAJECTORY AND ITS NEIGHBORHOOD UNDER 
CONSIDERATION (3 SUB INTERVAL CASE) 

For the piarpose of comparison of IDP and MLV methods 
for a one~hydro--one-thermal problem, the case in which the 
interval of optimization is subdivided into 3 equal sub- 
intervals is considered (see Figure 4.1). The points a, 
b, c, d represent the nominal state trajectory for a 
particular iteration number. The states represented by e, 
f, g, h indicate the neighbourhood considered in these 
methods. In the IDP method, 9 trajectories viz. aegd, aecd, 
aehd, abgd, abed, abhd, afgd, afed and afhd, need to be 
compared to find the optimal trajectory in this neighbour- 
hood assuming of course that all these trajectories are 
feasible. In MLV, 5 trajectories are to be compared in the 
worst case, and 3 trajectories in the best case. This 
result is obtained as follows: For the worst case: 

(l) abed is compared with aecd resulting in failure, that 
is resulting in the answer that abed is a better trajectory 
than aecd , (2) abed is compared with afed resulting in 
success ( even if this results in failure the ultimate 
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resuli: will be the saae) (3) afcd is compared with afgd 
resulting in failure ( 4 ) afcd is compared with afhd 
resulting in success (or failure}. Thus a maximtan nianber 
of 5 trajectories need to be compared in the MIV. In a 
Similar way it can be shown that in the best case only 
5 trajectories need to be compared in this method. 

In each iteration of IDP the neighbouring alternatives 
to the existing discharges are determined by measuiring the 
so called 'weighted output* of the hydro system [lO]* Thus 
in each iteration the quantized state values, the corres- 
ponding weights and the corresponding discharges 

need to be determined and stored to locate the best 
trajectory in the neighbourhood. Thus for the three sub- 
inteival problem, 24 main storage locations are required 
for IDP, In MLV, it can be seen that 15 storage locations 
are required (4 for the nominal states a, b, c, d, 3 for 
the nominal discharges, 3 for the incremental performance 
index values in the three subintervals, 1 for the perturbed 
state, 2 for the two discharges in the two adjacent sub- 
intervals of the perturbed state, 2 for the incremental 
performance index values in these two sub intervals corres- 
ponding to the perturbed state). 

In this comparison, the nmber of calculations^^^ ^^ ^ 
involved in either case is arrived at in the following 
way. In IDP while exploring all possible discharges which 
affect a transition from one stage to another stage, it is 
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observed ttLat only 5 distinct discharges are involved and 
need to he considered. Each of these distinct discharges 
results in a different weighted output of the hydro system; 
thus each distinct discharge is associated with a single 
calculation. In each of the first and last suhintervals, 
only two calculations are involved. Thus in all 9 cal- 
culations are involved per iteration in HP method. In 
the worst case of MLV if the variation from point h to e 
results in failure, the variation from point h to f is 
tried. Each perturbation involves 2 calculations, corres- 
ponding to the discharges in the two adjacent subintervals. 
Thus 4 calculations are involved in the perturbation pro- 
cedure referred to at the time instant k=2. A similar 
result holds for pertxir bat ions at the next time instant. 
Thus, in the worst case, 8 calculations are involved in 
MLV, Table 4.6 compares IDP and MLV for a one-hydro— one- 
thermal problem. It is assumed that identical procedures 
are used for calculations for the thermal subsystem in 
both cases. 

Prom Table 4.6 it can be seen that the MLV rec^tiires 
significantly less number of storage locations, less number 
of trajectories to be investigated and fewer number of 
calculations than IDP, 

The application of the MLV to higher dimensional 
problems,!,©. to problems of multireservoir hydro- 
■ thermal power systems, a procedure similar to the 
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Successive Approximation to D 3 mamic Programming (S,A.D.P.) 
[ 52 ] is adopted as can be seen in Sec, 4.7. In this method, 
the state trajectory of one hydro reservoir is subjected to 
the MLV, keeping the trajectories of all the other hydro 
reservoirs at their nominal values determined in the 
previous iteration. After completing the variation of one 
trajectory over the interval of optimization, the variation 
of the trajectory of the second hydro reseivoir is carried 
out in the same way. Thus the trajectories of all the 
hydro reservoirs are successively subjected to the MIV to 
complete one iteration. As it can be seen that an 
n-dimensional problem is split into n one-dimensional 
problems and hence the computational storage requirements, 
number of trajectories to be explored and number of cal- 
culations to be performed etc, vrill vary linearly with the 
number of the hydro stations imder consideration. Similarly 
in the case of IDP, adopting the S.A.D.P, procedure makes 
the requirements vary linearly with the number of hydro 
stations. Hence the comparison given for one-hydro-one- 
thermal problem given in Table 4*6 is also applicable to 
the multidimensional case. 

Thus it can be seen that for the same neighbourhood 
size, MLV is preferable to LDP, since the storage required 
and the calculations involved are less in the former case. 
When the initial trajectory is far away from the optimal 
one, ML? provides a faster transition from the current 
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nominal "trajectory to the next one compared to IDP, If . 
the optimal trajectory is encased in the neighbourhood at 
any stage of the iterative procedures of IDP and MI7, the 
IDP leads to the exact optimal trajectory, whereas M17 
provides a trajectory which yields a reduction in the 
value of the performance index compared to the one 
corresponding to the existing trajectory. Thus it can be 
seen that it is advantageous to use MLV initially, till 
we are close to the optimal trajectory and then switch 
over to IDP for convergence to the optimal trajectory. 
Alternately the difference between the two optima yielded 
by IDP and MIV can be reduced by choosing smaller varia- 
tional step size in MIV as the optimum is approached, 

4.7 FORMULATION OP THEEE-HIDRO-POUR-THEEMAl PROBLEM 

In this section, the mathematical formulation is 
given for a poxier system consisting of three hydro and 
four thermal generating stations. The hydro stations are 
assumed to be of the storage reservoir type, which are 
independent of each other. The hydro system is described 
by the state eq.uations 

x^(k+l) = ~ - ©2.(15:) (4.11) 

X2(k+l) = 2;2(k) + L2(h) - U2(k) - e2(k) (4.12)^^^^^^^^^^^^^^^^ 

x^(k+l) = x^(k) + L^(k) - u^(k) - e^(k) (4.13) 

In eqns. (4.1l) - ( 4.13) , h varies from 1 to N and the 
* initial and final state value x^(l ) and x^(N+l), respectively 
(i=sl,2,3) are specified. The performance index 
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■ N 4 

J - I I P.(P..(k)) (4,14) 

1&=1 j=l ^ 

is to "be minimized subjeoi to tlie ectnality coastmaint 


4 3 

PjjCk) = j|:i ^ "■ ^Loss^^^' (fc=l,...,H) 

(4.15) 

and the following inequality constraints 


min— -- max 

min ^ - =^1 ma.x 


(i=l, ...,4),(fc=l,,^,,,N) 

(4.16) 

( i=l, , . . , 3 ) , (lc=.2, , . . ■jIT) 

(4-.17) 


u > . < 

t min — 


u^(k) 


"'^i max ^ 


(i=l,,. .,3),(k=l,...,lT) 


(4.18) 


In eqn.(4.14) 


PjUjjCk)) =aj 

Similarly in eqn, (4.15) 

PjiCls:) =-^[1 + -^(zj^(k) + x^(k+l)) jUj^(k) 

^Losat'^) = X il 

In •eqn,(4.2l) 




P^(3s) ^ P^^(k) 


f ( i— ly • , . y 3 ) 

y ( 1 = 4 y * * . j T ) y ( 0 — 1 y • . . y 4) 


(4.19) 

(4.20) 
(4-21) 

(4.22) 


In the MlVy for the nominal trajectorjr at any stage 
of the iterative proced-ure,. the corresponding dis charge 
rates u^(k), and hydro generations ' Sox all k can 



•be calculated from eq.ns. (4.11) to (4,13) and (4.20) 
respectively. To determine the corresponding optimal 
thermal generations, it is necessary to solve the 
following mathematical programming problem in each 
sub interval represented by k. 


4 

Minimize I l^.(P^^.(k)) (4.22) 

subject to the equality constraint 


Pj5(k) 




d=i 


and the inequality constraints 


(4.23) 


min- ^Tj max ’ (d-lf •;4) (4.24) 

The Augmented Penalty Function Method, which has been 
discussed in detail in Chapter 3 has been chosen to 
solve the above problem. 


4. 8 ALGORITHM PBSORIPTION 

Step 1: Starting with a nominal trajectory for each 
hydro reservoir, the feasibility of the trajectory is 
established as described in Sec. 4.3. If the initial 
trajectory chosen is not a feasible one, then the 
algorithm automatically changes the trajectory to 
obtain a feasible one. 

Step 2: The discharge rates and hydro generations are 
cal.culated in each subinterval for the nominal trajec- 
tory. 
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Step 5: With the hydro generations Iniown over each 
of the sub intervals , a mathematical programming problem 
stated in eqns. (4.22) to (4.24) is solved repeatedly 
for all the subintervals to determine the corresponding 
optimal thermal generations. The subinterval costs, as 
well as the total cost over the entire interval of 
interest is also computed and stored. 

Step 4: Taking one hydro station at a time (in some 
order), the initial trajectory is varied by a small 
increment Ax, in either direction using the MI V at 
time instant k=2. feasibility of this variation is 
established as in step 1 above. Then the correspon- 
ding optimal thermal generations are determined in 
sub intervals 1 and 2 by solving the mathematical 
programming problem for these two sub intervals , The 
performance index value for these two subintervals 
is computed and compared with the corresponding value 
for the unvaried trajectory. If the variation consi- 
dered results in a reduction in the performance index 
value, the variation is a success and the process is 
shifted to the instant k=3. Otherwise the variation 
- A X is tried. If neither of these variations results 
in a reduction in the cost functional for the two sub- 
intervals indicated, the variational process is shifted 

to the instant k=3, retaining the original value of the 

■ » 

component at instant k=2. After the process is completed 
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at all instants k=2,...,IJ' for tlie first hydro trajectory, 
the procediire is snccessively repeated with the second and 
third hydro trajectories at time instants k=2,... ,]?. It 
is to he noted here that while the variation process is 
applied to one hydro trajectory, the other hydro trajec- 
tories are hold imchanged at their nominal values deter- 
mined in the preceding iteration. This completes one 
iteration. Thtis each iteration generates a feasible 
trajectory which is better than (or at least as good as) 
the trajectory in the previous iteration. 

Step 5 ° The above iterations are continued until a 
satisfactory convergence to the optimum cost is attained. 
This is judged from the difference in performance index 
values for two consecutive iteration numbers; this 
difference must be less than a prespecified small value 
for satisfactory convergence. 

4. 9 EXAMPLE 

Three— Hydro- Four -Thermal Problem; 

Results will be presented for the constant as well 
as variable step size as in Sec. 4.5. 

Hydro Reservoir 
_ ji” III” 

Maximirm storage (initial 

storage) in meter^ /sec.— month 80.0 135»0 80.0 

Final storage at the end of 

one year in meter ^/sec.— month 80.0 135 *0 80.0 



Hydro Reservoir 
I II III 


Minimum storage in meter^/ 


sec. -month 

0.0 

0.0 0.0 

Maximum allowable discharge 



rate in meter ^/sec.- 

73.0 

75.0 80.0 

Minimum allowable discharge 



rate in meter^/sec. 

3.334 

4.353 3.543 

Maximum hydro generation 
in MV 

89.03 

57.61 75.0 

Minimum hydro generation 



in MV 

O 

o 

0.0 0.0 

Basic head in meterSjH^j^ 

98,0 

50.0 75.0 

Constant c^^ 

0.004 

0. 002 0. 004 

Optimization period considered is one year of 12 monthly 

sub intervals . 



Constant G = 100.0 



Transmission loss Coefficients: 

= B 22 = 0.0005 



= B, . = 0.0008 

44 


®55 

= Bgg = 0.0007' 


®77 

= 0.0009 and 



= 0.0 (if^j) 

Thermal Station Characteristics: 



Thermal Station 


I 

II III IV 

a . in j8/MV-hr. 

2.5 

2.6 2.8 2.0 

b^. in ^/MV- hr. 

0.05 

0.06 0.08 0.08 

Maximum Thermal genera,tion 
in MV 

50.0 50.0 5S).0 50.0 

Minimum Thermal generation 
in M¥ 

0.0 

0.0 0.0 0.0 
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Hydro Reservoir 



I 

II 

III 

Minimum storage in meter^/ 




sec. -month 

0.0 

0.0 

0.0 

Maximum allowahle discharge 




rate in meter ^/sec.- 

75.0 

75.0 

80.0 

Minimum allowable discharge 




rate in meter^/sec. 

3.534 

4.353 

3.543 

Maximum hydro generation 
in M¥ 

■89.05 

57.61 

75 .0 

Minimum hydro generation 




in MW 

0.0 

0.0 

0.0 

Basic head in master s,Hq^ 

98.0 

50.0 

75.0 

Constant c^ 

0.004 

0.002 

0.004 


Optimization period considered is one year of 12 monthly 
suhintervals . 

Constant G = 100.0 

Transmission Loss Coefficients; = ^22 ~ 0.0005 

B,, = B. , = 0.0008 
00 44 

B 55 = Beg = 0.0007' 

Br^^ = 0.0009 and 

B^^. = 0.0 (i^j) 

Thermal Station Characteristics: 


a in ^/Mtf-hr. 

h^. in hr . 

Maximum Thermal generation 
in M¥ 

Minimum Thermal generation 
in M¥ 


Thermal 
I II 

2.5 2.6 

0.05 0.06 

50.0 50.0 

0.0 


station 


III 

IV 

2.8 

2.0 

0.08 

0.08 

55). 0 

50,0 

0.0 

0.0 


0.0 



TABLE 4.7; INELOW DATA FOR THE HYDRO RESERVOIRS 
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00 
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0 
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0 
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neglected in this prohlem. 
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COFSTAITT STEP SIZE CASE: ( Ax = 0.5 meter Vs ec. -month 
for all the hydro reservoirs). 

Table 4.8 shows the initial and optimal trajectories for 
the three hydro reservoirs. 


TABLE 4.8: INITIAL AND OPTIMAL TRAJECTORIES OF HYPRO 
RESERVOIR STORAOES 


Instant 

Ir 

St orage 

in the Hydro Reservoir in 

meter ^/s 

ec.-mon' 

ix 

Hydro 

reservoir 

I 

Hydro 

reservoir 

II 

Hydro reservoir 
III 


Initial 

Optimal 

Initial Optimal 

Initial 

Optima] 

1 

80.0 


80.0 

135.0 

135.0 

80.0 

80.0 ; 

2 

T5.0 


73.0 

128.0 

128.0 

75.0 

75.0 i 

3 

70.0 


70.5 

123.0 

123.0 

70,0 . 

.69.5 i 

4 

65.0 


64.5 

120.0 

119.5 

65.0 

65-0 1 

5 

60.0 


67.5 

117.0 

104-.5 

60.0 

66.5 ! 

6 

58.0 


65.5 

114.0 

99.0 

58.0 

64.0 i 

7 

53.0 


61.0 

110.0 

103.0 

53.0 

59.0 I 

j 

8 

59.0 


59.0 

115.0 

87.5 

59.0 

57.0 1 

9 

68.0 


70.0 

119.0 

89.5 

68.0 

69.0 

10 

72.0 


80.0 

122.0 

135.0 

72.0 

80.0 

11 

77.0 


80.0 

128.0 

135-0 

77.0 

80.0 

12 

78.0 


78.0 

131.0 

131.0 

78. 0 

78.0 

W 

80.0 


80. 0 

135.0 

135.0 

80.0 

80.0 


TaTpla. 4 . 9 shows the optimal discharges for the hydro reservoi 
over the entire scheduling interval. Table 4.10 gives the 



optimal schedules for hydro and thermal stations along wi 
the average load demand on the system. Tables 4.11 to 4. 
provide the optimal reservoir storage, discharge rates ar 
optimal schedule of hydro and thermal generations respec- 
tively using a variable step size procedure discussed in 
Sec. 4 . 4 . 

TABLE 4.9; OPTIMAL DISCHARGE RATES FOR HYDRO RESERVOIRS ; 


Sub interval 
k (month) 

Optimal 

discharge rate 

in meter^/sec. 

Hydro 

reservoir 

Hydro 

I reservoir II 

Hydro 

reservoir III 

* 

1 

7.8 

30.0 

5-8 

2 

12.5 

22.0 

15.5 

3 

17.0 

15.5 

15.5 

4 

8,9 

25.0 

10.4 

5 

19.0 

25.5 

19-5 

6 

23.0 

36.0 

23.5 

7 

29.6 

68.5 

29.6 

8 

32.8 

61.0 

31*8 

9 

46.4 

16.5 

45.4 

10 

40.3 

48.0 

40.3 

11 

32.1 

59.0 

32.1 

12 

44.3 

44.0 

44.3 
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VARIABLE STEP SIZE: ( Ax = 80.0 meter ^ /sec. -month for all 
the hydro reservoirs). 

TABLE 4.11: INITIAL AUD OPTIMAL TRAJECTORIES OP HTLRO 
RESERVOIR STORAGES 

Instant Storage in the Hydro Reservoir in meter ^/sec .-month 

k 

Hydro reservoir Hydro reservoir Hydro reservoir 
I II III 

Initial Optimal Initial Optimal Initial Optimal 


1 

80.0 

80.0 

135.0 

135.0 

80.0 

80.0 

2 

75.0 

65.0 

128.0 

123.0 

75.0 

65.0 

3 

70.0 

70.0 

123.0 

105.0 

70.0 

70. 0 

4 

65-. 0 

65.0 

120.0 

97.0 

65.0 

65.0 

5 

60.0 

67.5 

117.0 

96.5 

60.0 

62.5 

6 

58.0 

63.0 

114.0 

105.0 

58.0 

60.5 

7 

53.0 

55.5 

110.0 

85.0 

53.0 

58.0 

8 

59.0. 

56.5 

115.0 

87.12 

59.0 

56.5 

9 

68.0 

69.25 

119.0 

134.81 

68.0 

68.31 

10 

72.0 

79.81 

122.0 

134.87 

72.0 

79.81 

11 

77.0 

79.81 

128.0 

131.0 

77.0 

79.81 

12 

78.0 

78.0 

131.0 

131.0 

78.0 

78.0 

13 

80.0 

80.0 

135.0 

135.0 

80.0 

80.0 
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TABLE 4.12: OPTIMAL DISCBA.RGE RATES FOR HYDRO RESERVOIRS 


Sub interval 
k (month.) 

Optimal discharge rate 

in meter^/sec. 

Hydro 

reservoir I 

Hydro 

reservoir II 

Hydro 

reservoir III 

1 

15.8 

30.0 

15.8 

2 

5.0 

22.0 

5.0 

5 

16.0 

30.0 

16.0 

4 

9.4 

18.0 

14.4 

5 

21.5 

20.5 

19.0 

6 

26.0 

31.5 

21.0 

7 

26,6 

73.0 

29.1 

8 

51.0 

60.9 

32.0 

9 

45.8 

14.3 

44-. 9 

10 

40.3 

47.9 

40.3 

11 

31.9 

58.9 

31.9 

12 

44.3 

44.0 

44.3 





TABLE 4.15: OPTIML SCHEDULE FOR THE THREE-HTDRO-POUR-TIIERMAL PROBLEM 
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Table 4.14 shows the effect of variable step size on the 
computational time. 

TABLE 4.14: EPPECT OP THE STEP SIZE Ax OH THE GOMPUTATIONAI 
TIME 


Step size Ax in 
me ter^/sec. -month 
(for all the three 
hydro stations) 

Cost in 

Initial 

dollars 

Pinal 

Computational time in 
seconds (Execution) 

80. O'*" 

7082.77 

6976.97 

165 

20.0+ 

7082.77 

6976.97 

165 

1.0 

7082.77 

6992,75 

205 

0-5 

7082.77 

6991.77 

208 


Note: + in the above table indicates a variable step size'. 

4.10DIS0USSI0ir 

Numerical results are presented in Tables 4.1 to 
4*5 for a one-hydro-one-thermal system and in Tables 4.7 
to 4.14 for a three-hydro-four-thermal system respectively. 
He suits are presented for both the constant step size and 
variable step size separately in both the cases. It can 
be seen from the numerical results that the variable step 
size procedure resulted in a faster convergence as 
compared to the constant step size case. The initial 
value in the variable step size procedure- can be as large 
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as possilDle (note that it can infact he equal to the "value 
of the maximuia allowable storage in the reservoir itself;, 
however, in a mult ires ervoir system, it implies a value 
equal to the mazimum allowable storage in the smallest 
reservoir)* 

fables 4*5 and 4.14 summarize the effect of the 
variable step size on the computational time consumed in 
both the problems respectively. The computational time is 
significantly less with variable step sizes compared to the 
constant step size. Thus on the basis of the experience 
gained, it is recommended that the initial step size in 
MIV procedure should be as large as possible and then 
progressively reduced as the optimum is approached. 

4.11 CONCIUSIOFS 

In this chapter, the long range optimal power 
generation schedule in a hydro thermal power system is 
obtained using an approximate model discussed in Sec. 2. 2. 
Method of Local Variations has been employed to modify . 
the hydro trajectory in each iteration, starting from an 
initially assumed trajectory to reach the optimum. 

Augmented Penalty Function Method has been used to solve 
the mathematical programming problem encountered in each 
subinterval for solving the optimal thermal generations. 

The Method of Local Variations has been found to 
be attractive due to its simplicity, efficiency and ease 
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of implementation, limits on the. hydro reservoir, the 
hydro and thermal generations have "been efficiently taken 
care of by this method. Furthermore, the method requires 
less storage and seems to require less computational time 
than those for the Incremental Dynamic Programming 
procedure adopted by Bernholtz and Graham to solve the 
same optimization problem. A recommendation based on 
computational experience is made regarding the choice of 
the step size while working witii the MLY algorithm. For 
faster convergence, it is recoliimended that a step size as 
large as is consistent with the problem is selected to 
start with and gradually reduced as the optimiun is 
approached. This provides a significant reduction in 
computational time compared to the use of constant step 
size throughout the process. 



CHAPTER 5 


SOLDTIOH TO THE SHORT RANGE SCHEDULING PROBLEM 

5.1 INTRODUCTION 

In Chapter 4, the problem of long range optimal 
scheduling in hydro-thermal power systems is solved 
using an approximate model, wherein the characteristics 
of the eleotrical network have been represented by the 
loss coefficients and the limits on the voltages, re- 
active powers, line flows etc. are not considered. 
However, in a short range scheduling problem when the 
data available would be more precise, in addition to 
hydro system dynamics one would like to represent the 
electrical network more rigorously using the power 
flow equations, and consider the various constraints 
on voltages, reactive powers, line flows etc. Such a 
sophisticated model has already been described in 
Sec. 2.3. Thus the attention, in this chapter is 
focussed to the solution of the problem using this 
model. The head variations due to the reservoir 
dynamics have also been taken into account for the 
sake of generality. 

To the author’s knowledge, so far only two 
attempts [35], [20] have been made in the past to 
solve the short range scheduling problem both using 
A.G. power flow model for the electrical power network 



as in Sec. 2.5. Howe-ver the- hydro system dynami-os has 
not been included in [35]. Both the above attempts are 
briefly summarized below. Ramamoorty and Gopala Rao 
[53] have formulated this problem as an additively 
separable nonlinear programming problem. The repre- 
sentation of the hydro system is as follows. It is 
assmed that the total energy available at each hydro 
station is known for the total interval of optimization. 
Under the assumption of constant head and equal length 
of each subinterval, it was possible to write a coupling 
constraint for the total interval of optimization for 
each hydro station, which stipulates that the available 
energy at each hydro station is utilized in power 
generation. Associating a dual variable (which has the 
physical significance of the equivalent water cost) 
with the coupling constraint at each hydro station, the 
problem thus reduces to minimizing an objective function 
(being the sum of the cost of thermal generation in all 
the' subintervals ) , coupling equations defined by the 
total energy constraints as seen above, and the 
constraint sets defined by the load demand and the 
various operating constraints on the system in the 
corresponding subintervals. With a choice of the 
above mentioned dual variables the problem of IT sub- 
intervals was split into N subproblems, one for each 
subinterval (which can be solved separately) ■using 
lasdon's [22] decomposition technique. The solution 
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procedure consists of solving these suhprohlems 
iteratively using S'®'!! [35], while updating the dual 
variables in each iteration by a gradient technique 
until the convergence to the optimum solution is 
obtained. 

As can be seen from the above discussion, no 
attem,pt was made by them to represent the dynamics of 
the hydro subsystem in their formulation, Thiis there 
was no method of checking the operational constraints 
of the hydro subssrstem, such as the upper and lower 
limits on the hydro reservoir storages and the associa 
ted discharge rates etc. Also the relaxation of the 
assumptions of constant head and equal length of each 
subinterval will make it very difficult to write the 
above mentioned coupling constraints for the hydro 
stations. Furthermore, unless the initial choice of 
the dual variables is proper, convergence to the 
optimum may be difficult in addition to the large 
computation time required. In the author’s opinion, 
the hydro subsystem was not represented adequately 
.and their method of solution was the same as that of 
solving the optimal scheduling problem of a purely 
thermal system except that the additional term in the 
cost function representing the equivalent cost of 
hydro power. 



82 


Bonaert, El-ABiad and Koivo [20] have solved 
the short range schediiling problem using a model 
similar to the one considered in this chapter. Their 
approach was basically the same as in [10]. The 
solution procedure in each iteration consists of 
essentially starting with a nominal trajectory of 
reservoir storages and applying the Dynamic Programming 
algorithm locally in a neighbourhood around the nominal 
trajectory and in the process to obtain the optimum 
trajectory. In each sub interval, knowing the hydro 
generations from the hydro subsystem, the corresponding 
optimal thermal generations are obtained by solving an 
optimal power flow problem. At the end of each itera-. 
tion the existing nominal trajectory is replaced by the 
best trajectory in its neighbourhood. The difficulty 
with this method as seen from the numerical results 
reported [20] was that the computational storage re- 
quirements were significantly high even for the 5 bus 
system considered by them. 

In this chapter, the solution to the above 
problem is obtained employing a similar decomposition 
of the combined system into hydro and thermal subsystems 
as that of Bonaert etal. The major difference between 
these two approaches stems from the method of modifying 
the initially chosen nominal trajectory in successive 
iterations. In the present method the MLV which has 
been applied in the long range sched\xling problem in 
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Chapter 4 is also used for the short range scheduling 
problem. In each iteration, starting -with a nominal 
trajectory of reservoir storages, the method consists 
of perturbing the initial trajectory at one instant at 
a time and choosing the best trajectory in two adjacent 
sub intervals comparing the three possible alternatives. 
Dtiring the perturbation at any instant, the rest of the 
trajectory is assumed to remain fixed at its previous 
value. The perturbation process is performed at all 
instants, to complete one iteration. As in [20], 
knowing the hydro generations from the hydro subsystem • 
in each subinterval, starting from the end of first 
subinterval to the beginning of the last subinterval, 
the corresponding optimal thermal generations are 
calculated. In Sec. 4.6, it has been shown that the 
storage regttiirements for the MI¥ are significantly less 
than those for the IDP. Thus it can be seen that the 
present method affords a better approach than that of 
Bonaert etal from the aspect of computer storage 
requirements. It can also be seen from the numerical 
results reported later in this chapter, that the 
computational time requirement is also less compared to 
that of Bonaert etal. 

5.2 PROBIBM BEPIlSriTIOH 

The system configuration considered for short 
range scheduling problem in this chapter is as shown 
in Pigure 5.1. This is similar to that discus sed by 



Hydro station 2 


BA 



Hydro station .1 

FIGURE 5.1: HISTSM S'TUEIED FOR THE SHORT RAHGE RROBIEM. 


Bonaert etal [20] and the motivation for this is so 
that a comparison between the present approach and that 
of Bonaert etal can be made. Kie generation at buses 
1 and 4 are from hydro stations, which are assumed to 
be operating independent of each other and thermal 
generations are assumed at buses 2 and 5. loads 
It, 1- and L„ are assumed to be at buses 5, 2 and 5 
respectively. The detailed mathematical model for a 
general short range problem has been described in 
Sec. 2.3. let the total interval of optimization 
(one day) be divided into 6 equal subintervals of 
4 hours each duration. It is assumed here that the 
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loads in each, suh interval remain constant at all 
load buses. The initial and final storages of the 
reservoirs are specified. 

The state equations describing the hydro 
subsystem dynamics are 

x^(k+l) = x^(k) + L^(k) - u^(k), (i=l,2), (k=l,...,6) 

. (5.1) 

x^(l) and x^(6), (i=l,2) are specified. The hydro 
generation at the i-th station during the k-th sub- 
interval can be related to the water head and discharge 
rate by the relation [16] 

Hnn C-; 

= -^[l + -^(x^(k) + Xj^(k+l)) ]u^(k) , 

(i=l,2), (k=l,...,6) (5.2) 

Bounds on the reservoir storages and discharge rates 
are 

"i min < i mam ’ (k=l,...,6) 

(5.3) 

min - max ’ (^ = l>2)i (k=2,...,6) 

(5.4) 

The thermal subsystem can be formulated as 
follows (from Sec. 2.5). In each subinterval, having 
fixed the hydro generations at the two hydro stations 
from the' nominal trajectory, it remains to solve ah 
optimal power flow for each subinterval. Hence this 
problem becomes a static optimization prob3.em in each 
sub interval which can be stated as 
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6 2 

Minimize J = I J F (I>rj,.(k)) (5.5) 

k=l i=l 

where J represents the total thermal generation cost 
subject to the following equality and inequality 
constraints . 


Equality Constraints: 

Ihe real power balance is given by 

Ii(k) - P^(k) + C^(k) = 0 , (i=l,2), (k=l,....,6): 

(5.6) 

The reactive power balance is given by 


K^(k) - Q^(k) + D^(k) = 0 , (i=l,2), (k=l»...,6) 

(5.7) 

where the real and reactive power injections at the 
i-th bus are respectively 


Ii(k) 


V^(k)?^(k)T^^Sos(0. .+6^(k)-6^(k)) 


+ v2(k)Y^^eos9ij. , (k=l,...,6) (5.8) 


Ki(k) 


I F^(k)Vj(k)I.jSin{6ij+6i(k)-6j(k)) 

+ ^i(l!:)YiiSinei^ , (k=l,...,6) (5,9) 


Inequality Constraints : 


The various operating restrictions in the system 
impose the following limits on the variables P^(k), 
Qj_(k) and Vj_(k). • 
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o 

Vl 

(5.10) 

min - 

< 0 

(5..11) 

^i min ** 

< 0 

(5.12) 

- iq 

'< 0 

(5.13) 

h min - 

< 0 

(5.14) 

\ max 

£ 0 

(5.15) 


The restrictions on the individual line flows to 
ensure system stability can be represented by 

|6i(k) - 6j(k)| < |T,3 (5.16) 

Eqns. (5.10) to (5.16) are to be satisfied in each 
sub interval. Thus solving the thermal subsystem for 
the optimal thermal generations for known hydro 
generations consists of minimizing (5.5) subject to 
equality constraints (5.6) and (5.7) and the upper 
and lower limits on the problem variables given by 
(5.10) to (5.16). 


5.3 TWO IE7EL SOLUTION APPROACH 


The combined system is decomposed into a 
dynamic hydro subsystem and a thermal subsystem { 
(static in each subinterval) , In the hydro subsystem, 
starting with a nominal (feasible) trajectory of 
reservoir storages, the nominal hydro generations in 
the various sub intervals are obtained. Coiasidering 
these hydro powers as injections into the electrical 
network at the appropriate buses, the thermal subsystem 
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in each subinterval is optimized using the Optimal Power 
Plow Solution due to Pommel and Tinney discussed in 
Sec, 3.4. The existing nominal hydro trajectory is now 
varied using the MLY iteratively to result in a reduc- 
tion in the cost of thermal generation in each iteration, 
while an optimal power flow is obtained with the given 
hydro injections, for the thermal subsystem, in each 
subinterval. As can be seen, this solution procedure 
constitutes a two level approach, to coordinate the 
hydro power iteratively to result in an optimal schedule 
for the combined system, while the two subsystems are 
decoupled at the first level with the hydro power 
injections as the interaction variables. Figure 5.2 
shows the flow chart of .the algorithm, which is 
described in detail in the next section. 

5.4 ALGORITHM DESCRIPTIOH 

Step 1; To start with an initial nominal trajectory 
(k=l,...,6), for all hydro stations, (i=l, 2) 
is chosen. Feasibility of the trajectory is established 
by making sure that the trajectory does not violate the 
constraints on states x^(k), and the resulting discharges 
Ujj_(k), for all i and k. If the initial trajectory chosen 
is infeasible, the program changes the trajectory to 
obtain a feasible one as described in Sec. 3.2. Knowing 
the inflow data for all the reservoirs over the interval 
of optimization, the sub interval nominal discharge rates 
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FIGinffi 5.2? OPTIMAL SCHEDDLIHff ALGORITHM POS 
THE COMB UTEL STSTM . 
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•unchanged, while the trajectory of one of the hy<3.^^ 
station is "varied, Peasibility of this variation 
established as in step 1 above. The discharge 
and the hydro generations for subintervals 1 and ^ 
are obtained using eqns. (5.1) and (5.2) and the 
corresponding optimal thermal generations are det®^*” 
mined by obtaining the optimal power flow solutiof*-^ 
for these two subintervals. The subinterval costs 
for these two subintervals are obtained and 
with the corresponding two subinterval costs for hli® 
unvaried trajectory. If the variation considered 
results in a reduction in the cost functional for 
two subintervals indicated, the variation is a succ^^^ 
and the process is shifted to the instant k = 3. 
Otherwise the variation - is tried. If neither 
these variations resul'fcs in a reduction in the cost 
functional for the two subintervals considered, the 
nominal trajectory value is retained at this time 
instsnt and the procedure is repeated at time 
instants k = 3,4,,.. etc. When this procedure has 
been applied to the states at all the time instants 
k = 2,..., 6, the procedure is repeated with the 
second hydro station. When both the hydro stations 
have been tried, one iteration is said to have been 
compl6t6d , 

Step 4: The above iterations are continued tint il a 
satisfactory convergence to the optim-um cost is ^ ^ 
attained. 
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5.5 NIMERICAL EMPIB 

The algorithm described in the above section 
uses throughout a constant step size vector for 
modifying the hjrdro trajectory in each iteration. 
Since, it has already been established in Chapter 4, 
the-t starting with as large a step size vector as 
possible initially and then gradually reducing it 
during the iterative process as the optimum is 
approached, results in a significant reduction in 
the computational time, in this section, nmerical 
results obtained by variable step size vector alone 
are reported. 

The data as taken from [20] is as follows. 

The system considered is composed of the standard 
five bus system proposed by IEEE (Figure 5.1). The 
line data and the load data are given in Tables 5.1 
and 5.2 respectively. The data pertaining to the 
hydro and thermal systems are given in Tables 5.3 and 
5.4. Each sub interval considered is of 4 hours 
duration. Table 5.5 gives the upper and lower limits 
on the reservoir storages at the various instants over 
the optimization interval. Table 5.6 stipulates the 
limits on the voltages, active and reactive powers at 
the various buses. Table 5.7 gives the initial and 
optimal trajectories of hydro reservoir storages. 

Table 5.8 gives the optimal schedule, for hydro and 
thermal stations obtained using an initial step size 
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of A 2 : = 160.0 meter^/sec.-4 hours for "both the 
hydro stations and gradually reducing it as described 
in Sec. 4.4. Base power is taken to be 100 MVA» Bus 
number 5 is taken as the slack bus. 



TABLE 5.1: LINE 

LATA (P 

.N.) 


Line 

1-5 

2-5 

3-5 

2-3 

3-4 

Resistance 0.03 

0.08 

0.105 

0.033 

0.106 

Reactance 

0.103 

0.262 

0.347 

0.118 

0.403 


TABLE 5.2: 

ACT I'VE 

LOAN DATA (P.U. ) 


interval 

1 2 

3 

4 

5 

6 

load Lq^ 

0.6 0.8 

1.0 

0.9 

0.7 

0.5 

Load L 2 

0.3 0.4 

0.5 

0.45 

0.35 

0.25 

Load L^ 

0.3 0.4 

0.5 

0.45 

0.35 

0.25 


Note; Reactive loads are 10 percent of active loads. 



TABIE 5.3: CHARACTERISTICS OP THE HYDRO SUBSYSTEM 


Hydro system 

Hydro 1 

Hydro 

2 


(i=l) 

(i=2) 

Initial storage in 
meter^/sec,-4 Hours 

2330 

2330 


Pinal storage in 
meter^/sec.~4 Hours 

2350 

2330 


Maximum allowaHle discHarge 
in meter^/sec. 

500 

500 


Minimum allowaHle discHarge 
in meter^/sec. 

100 

100 


Initial Head of tHe reservoir 



in meters 

14 

14 


Pinal Head of tHe reservoir 




in meters 

14 

14 


Basic Head H^^^^ 

10 

10 


Inflow (assumed constant 
tHrougHout in meter^/sec. 

250 

300 


TABLE 5.4: THERMAL SUBSYSTEM 

DATA (Same 

for HotH 

stations) 




Lower limit of power (p.u. ) 



0.1 

Upper limit of power (p.u.) 



1.0 

Linear coefficient of tHermal cost 



(^/M¥-Hr.) 



3.8 

Quadratic coefficient of tHermal cost 



(^/))I¥^-Hr.) 



0.01 



TABLE 5.5; UPPER ATTB LOWER LIMITS OR TilE RESBRYOIR STORAGES 




TABLE 5.6: EMITS ON VOLTAGES AlLD GENERATIONS 
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TABLE 5.8; OPTBIAL SCHEDULE EOR HIDRO AUD THERMAL STATIOITS 


Sub- 

int erval 
k 

1 

2 

5 

4 

5 

6 


0.2207 

0.2315 

0.2568 

0.3667 

0.3775 

0.3734 

PH2(k) 

0.3234 

0.3218 

0.3547 

0.3955 

O. 4 O 68 

0.4033 


0.4208 

0.6819 

0.8654 

0.6114 

0.3408 

0.1485 


0.2503 

0.3806 

0.5606 

0.4503 

0.3000 

0.1001 

u^(k) 

182.0 

190.0 

210.0 

300.0 

310.0 

3 O 8 .O 

U2(k) 

267.0 

265.0 

275.0 

325.0 

335.0 

333.0 


Note; All units are as defined before except u^(k) 
and U 2 (k) which are expressed in meter^/sec. 
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TaTole 5.9 gives the voltages, phase angles and reactive 
powers at all the huses over the entire interval of 
optimization (for the optimum condition). 

TABLE 5.9: VOLTAGES, PHASE AHGLES AND REACTIVE POWERS AT 


THE VARIOUS BUSES 


Sub- 

interval 

k 

1 

2 

5 

4 

5 

6 

h 

1.02 

1.02-' 

1.02 

1.02 

1.02 

1.02 

h 

1.679 

1.744 

1.899 

2.567 

2.632 

2.607 

h 

1.04 

1.04 

1.04 

1.04 

1.04 

1.04 


-0.368 

-0.831 

-0.690 

-0.495 

-0.167 

-0.554 


1.0271 

1.0245 

1.0214 

1.0225 

1.0249 

1.0277 

s 

0.029 

-0.80 

-1.111 

-0.422 

0.354 

0.485 


1.02 

1.02 

1.02 

1.02 

1.02 

1.02 


7.862 

7.001 

6.972 

9.034 

10.042 

9.852 

^5 

1.04 

1.04 

1.04 

1.04 

1.04 

1.04 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


-0.2581 

-0.2609 

-0,2674 

-0.2949 

-0.2975 

-0.2965 

Qj 

0.1587 

0.1854 

0.1983 

0.2020 

0.1827 

0.1777 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


-0.0808 

-0.0758 

-0.0682 

-0.0768 

-0. 0842 

-0.0897 


0.3584 

0.3689 

0.4035 

, / 

0.4400 

0.4558 

0.3998 


Rote: Q^,...,Q^ are the reactive powers at buses I,..., -5 

respectively. All units are in p.u. except the phase 

angles (6 's) , which are given in degrees. 
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Initial cost in dollars = 7805 « 76 
Pinal cost in dollars = 7779*36 

Table 5»10 shows the effect of variable step size on the 
computational time» In addition to the variable step size 
case presented in detail, three more cases (one variable 
step size and t^o constant step sizes) are summarized in 
this table. Table 5*11 draws a comparison between the 
proposed method and that of Bonaert etal from the aspect 
of core storage and computational time req,uirements. 

TABLE 5.10: BPPECT OP VARIABLE STEP SIZE OR 
COMPUTATIONAL TIME 


Step size ax in 
meterVsec.-4 hrs. 
(for both the 
hydro reservoirs) 

Cost in dollars 

Comput at ional 
time in secs, 
(Execution) 

Initial 

Pinal 

160+ 

7805.76 

7779.36 

177 

40+ 

7805.76 

7779*36 

177 

10 

7805*76 

7779.88 

265 

5 

7805.76 

7779.36 

513 


Note; tindicates a variable step size 
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TABLE 5.11: COMPABISQB OF THE COMPUTATIONAL TIME AND 

STOEAGE EEQUIEEMBNTS FOR THE PROPOSED METHOD 


WITH THAT 

OF BONAERT ETAL. 


Method 

Core storage 

Computational 


required 

time in secs. 


in word memory 

(Execution) 

Bonaert etals’ 

method (CDC 6500 
computer) 

47. 2K 

200 

Proposed method 

(IBM 7044 computer) 

17.29K 

177 


5.6 DISCUSSION 

As seen earlier in tliis chapter, both the present 
approach and that of Bonaert etal [20] use the Optimal 
Power Plow Method in solving the thermal subsystem for 
optimal thermal generations. Though the actual steps in 
the solution are different in both the approaches, it can 
be seen easily that the computer storage required for 
thermal subsystem in either case remains same. The 
storage primarily consists of the sub interval costs, total 
cost and the subinterval optimal thermal generations etc. 
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But in the hydro subsystem Bonaert etal have used 
the IBP and the MLV is used in this approach. As already 
seen, the basic advantage of employing the MLV approach 
is for the hydro subsystem calculations. In Sec. 4.6, 
as also in [ 56 ] both the MLV and LDP algorithms have been 
directly compared from the aspects of computational 
requirements and it was established that in modifying the 
hydro trajectory in each iteration, the MLV requires 
remarkably less computational storage compared to the 
IDP procedure. This in the author's opinion is the 
main advantage of the present approach over that of 
Bonaert etal. The computer storage and time requirements 
of both the methods are given in Table 5.11. It can be 
easily seen from this table that the computer storage 
requirements for this method are less than half that of 
Bonaert etal. Regarding the computer time, the calcula- 
tions are carried on CDC 6500 computer by Bonaert etal, 
whereas IBM 7044 computer was used in this case, which is 
a much slower machine. The time required to obtain the 
solution to the numerical problem in Sec. 5.5 is 177 secs, 
on IBM 7044 whereas on CDC 6500, the time of computation 
as given in [20] was 200 secs. This clearly shows that 
there is a ..siguificant saving in computational time also , 
with the present method. Both the saving in storage and 
computation time make this method more attractive especially 
for a large interconnected power system. 
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5.7 OONCIUSIOUS 

In this chapter? the short range optimal scheduling 
problem using a sophisticated model for a hydro thermal 
power system has been solved using the Method of local 
Variations and Optimal Power Plow Method, Ihe combined 
system h8,s been decomposed into a hydro subsystem and a 
thermal subsystem and a two level solution approach is 
employed. The hydro subsystem has been represented by 
its dynamics considering the deterministically known 
inflows into the reservoirs and the load demand. The 
thermal subsystem has been represented by the A.C. power 
flow model. This would enable the limits on the voltage 
magnitudes, reactive powers and the line flows etc. in the 
electrical network to be considered. The modification 
of the hydro trajectory is achieved in each iteration 
using the MIV, while the optimal thermal generations are 
obtained in each subinterval by the Optimal Power Plow 
Method. 

Results are presented for a 5 bus system with fwo 
hydro and two thermal stations. A comparative discussion 
between this method and that of Bonaert etal is presented. 
The requirements of core storage and computation time 
for this method have been observed to be significantly 
less than those of Bonaert etal for the same 5 bus system. 
Hence the proposed method would be more attractive for 
large interconnected power systems. 



CHAPTER 6 


OPTIMAL SCHBDULIHO OP PUMPED STORAGE PLANTS 

6.1 INTRODUCTION 

In tlie problems considered, in Chapters 4 and 5, 
it has been implicitly ass-omed that the hydro plants are 
operated purely as generating stations and no pumping 
operation has heen considered at these stations. The main 
function of a pumping cum generating (pumped storage) 
plant is to pump water from the back-hay reservoir to the 
fore-bay reservoir during the periods of low demand and 
use the water thus stored for generation purposes during 
the periods of high demand. To utilizie the capacity of a 
pumped storage plant in a combined hydro-thermal power 
system in the most efficient manner, it is necessary to 
perform an optimization study for the system, aimed at the 
determination of an optimal pumping and generating schedule 
for the pumped storage plant. The schedule should be such 
that the total fuel cost incurred in thermal generation 
during a specified optimization interval is minimum. The 
resulting schedule should of course conform to the system 
operatir^ constraints such as the limits on the storages 
of the reservoirs, discharge rates (pumping as well 
generating), the thermal generations, the voltage magnitudes 
at various buses, the reactive powers and the line flows 
etc. as already discussed in Chapters 2 and 5* 
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The problem of pumped storage optimal scheduling is 
different from the purely hydro generation scheduling 
problem due to the fact that in a pumped storage plant, the 
plant operates as a generating station in certain sub- 
intervals and as a pumping station consuming power in the 
remaining subintervals of the optimization interval. 

Yarious attempts [57] to [41] have been made in the past 
to study scheduling with pumped storage plants, however 
most of them have dealt with specific cases rather than 
problem of a general nature discussed in this chapter. 
Bernholtz etal [37] have given a method using the modified 
dynamic programming procedure [42], [43] for the optimum 
scheduling of Ontario Hydro’s Sir Adam Beck-Niagara 
generating station. For the problem considered, water 
from upper Niagara river enters the Beck fore-bay which is 
connected by a pumping cum generating station to the 
storage reservoir and to the tailrace by means of a 
generating station respectively. The discharge cum 
pumping rates of the pumping cum generating station as well 
as the discharge rates of the Beck generating station 
depend on the volumes of water available in both the 
storage reservoir and the Beck fore-bay. The discrete 
scheduling problem consists of determining N pairs of 
discharges through pumping cum generating and the 
generating stations over the N subintervals ofoptimi- 
zation leading from a specified initial state 
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to a specified final state ^(K+l)^ where E^, 

and Ep, refer to the volixmes of the water in the fore-hay 
and the storage reseirwoirs respectively, while maximizing 
a snitahly defined objective function satisfying the 
various constraints on storages and discharge rates etc. 

The solution procedure consists of applying the Dynamic 
Programming starting backwards from the last subinterval. 

At the beginning of each subinterval, a two way table is 
constructed with both the volumes in the fore-bay and the 
storage reservoirs as the variables. Each combination of 
these two volumes is checked for feasibility and is denoted 
by a feasible cell, if the constraints are all satisfied 
and there is at least one pair of feasible discharges to 
reach a feasible cell at the next instant* Prom every 
feasible cell at the beginning of the F-th sub interval to 
reach the specified final state, the possible discharge 
rates and the corresponding weighted outputs are calculated 
and recorded. This results in a table for the F-th instant 
giving the feasible cells with relevant information recorded 
in them. Similarly tables are constructed for all instants 
i = 1 ,...,F, working from F-th instant backwards. Prom the 
tables so constructed, the optimum schedule can be obtained 
by identifying the best possible route. The procedure, as 
can be seen, becomes quite unwieldy if there is more than 
one pumped storage station in the power system and the 
computational storage and time requirements will also 
become quite high. Bernard etal [58] have considered one 



107 


pumped storage plant in conjunction with a thermal plant 
and suggested an iterative procedure of obtaining the 
optimum pumping generating schedule starting with a nominal 
schedule. Their procedure is based on the following two 
results obtained using the calculus of variations; 

(i) Optimum utilization of the available water is 
obtained when the hydro generation is scheduled for 
specified subintervals such that the incremental fuel 
savings are equal for each corresponding subintervals, 

(ii) Additional economy may be obtained as long as the 
incremental fuel saving exceeds the incremental pumping 
cost. To start with, incremental pumping cost curves as 
well as incremental fuel saving cost curves are determined 
for each of the pumping as well as the generating sub- 
intervals (prespecified) respectively as a fuaction of 
volumes of water. Starting from a nominal schedule of 
pumping and generating, the various incremental fuel 
saving costs and the incremental pumping costs for the 
corresponding subintervals are obtained from the above 
curves. The method consists of systematically improving 
the hydro schedule ensuring the adherence to the above 
two results. The iterative process converges to the 
optimum schedule when the incremental fuel savings cost 

is equal to the incremental pumping cost for each sub- 
interval in the pumping phase. This proced we becomes 
difficult for a problem consisting of more than one 
pumped storage and one thermal plants, as the fuel 
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saving cost and. pumping cost curves for the corresponding 
suhintervals as a function of water volumes cannot he 
independently obtained as in the simple case seen above - 
Further in their model, the transmission losses in the 
system are not considered. 

Brainbridge etal [39] used a gradient procedure 
on the hydro system of the problem and the d:ynamic program- 
ming procedure combined with the use of coordination 
equations for the thermal system. The basic procedure in 
this method is to start with a trial schedule, and by an 
iterative adjustment of the outputs of the hydro plants 
in the direction of steepest descent, the corresponding 
outputs of the thermal plants are obtained using the 
combination of dynamic programming procedure and the 
solution of the linear coordination equations. The 
iterative process is repeated until a minimum cost is 
attained. Cobian [ 40 ] considered the problem of a pumped 
storage plant with three interconnected power systems. 

The problem is to determine the optimal power inter- 
changes between the -different power systems and the pumped 
storage plant. Each power system is represented by its 
cost characteristic and the pumped storage plant by its 
reservoir dynamics. Transmission losses are neglected 
in the formulation (however, he has suggested an approxi- 
mate way of accounting for the transmission losses in the 
interconnected power system). Dynamic Programming procedure 
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is used to oltain an optimum schedule. In each suhinterval 
of optimization, a number of mathematical programming 
problems are to be solved and the modified version of 
Rosen’s [44] gradient projection technique is employed to 
solve these mathematical programming problems. As admitted 
by the author in his paper [40], this procedure is not 
applicable to systems having more than two pumped storage 
plants due to excessive fast access computer memory require- 
ments. Rees and Larson [41] have also considered a pumped 
storage plant in combination with an interconnected power 
system. The diversity exchange contracts also are consi- 
dered in the formulation of the problem and the transmission 
losses are neglected. Successive Approximations to Dynamic 
Programming (S.A.D.P.), [ 52 ] is employed to solve the above 
problem. 

In all the above attempts, it can be seen that 
basically Dynamic Programming with some modifications in 
certain cases has been employed in solving the pumped 
storage scheduling problem under consideration. Because 
of the very nature of the solution procedure requiring 
excessive core storage and computational time, the attempts 
have been limited to small systems. In this chapter, the 
Methods of local Variations and the Optimal Power Plow 
used in Chapter 5 have been employed again to solve the 
problem under consideration. The superiority' 'of' the 
basic approach employed here regarding the computational^^^ 
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time and core storage requirements has already heen shown 
in Chapters 4 and 5. The inclusion of the pumped storage 
plants does not alter the claims made preiriously regarding 
the computational core storage and time requirements. 

6.2 PROBIBH POEMIIATIOR 

The optimal scheduling problem in a power system 
consisting of pumped storage plants in addition to 
conventional hydro plants and thermal plants is formulated 
as a discrete time optimal control problem as follows. Let 
pfh and r represent the number of pumped storage, conven- 
tional hydro and the thermal plants respectively. Let 
x:^(k) be the storage in the reservoir of the i-th hydro 
plant (in the case of a pumped storage plant, it refers 
to the storage of the fore-bay reservoir) at the beginning 
of the k-th subintei^al. Similarly let u^(k) be the dis- 
charge rate from the reservoir of the i-th hydro plant 
(in the case of a pumped storage plant, this refers to the 
pumping rate of water from the back-bay reservoir into the 
fore-bay reservoir). The total interval of optimization 
is assumed to be divided into R equal subintervals of unit 
length for simplicity (as in Chapters 4 and 5)* The initial 
and final storages of the reservoirs in the p-umped storage 
as well as conventional hydro plants are ass-umed to be 
specified. The inflows into the reservoirs of the pumped 
storage plants are neglected in this formulation without 
loss of generality. The hydro subsystem dynamics is 
defined by the equat ions 
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x^(k+l) = x^(k) - u^(k) , (i=l,...,p), (k=l,.,.,N') (6.1) 

X, (k+1) = X. (k) +l.(k) *-u.(k), (i=p+l, . . , ,p+li) , 

*1 

(k=l,...,N) (6.2) 

Xj_(l) and Xj^(F+l), (i = l,...,p+h) are specified. It may 
be noted here that u^(k) is negative in the pumping mode 
of the plant and positive in the generating mode. The 
pumping and generating powers can be approximately related 
to the water head, discharge rate (or pumping rate) by the 
following relations. 

Power required during the pumping operation, by the i-th 
pumped storage plant in the k-th sub interval is 

Pjj^(k) = + ^(x^(k) + x^(k+l))]u^(k) (6.3) 

and power generated by the i~th conventional hydro or 
pumped storage plant during the generation operation in 
the k-th sub interval is 

Pg^(k) = + •^(x^(k) + Xj, (k+l))]u^(k) (6.4) 

where - 11 ^ and are the efficiencies of the plant in the 

pumping and generating modes respectively. It may be 
noted here that in eqn.(6.3)» if ii^(k) is negative (for 
pumping operation), then power Pjj^(k), will also be 
negative. The discharge rates, pumping rates and the 
storages of the hydro reservoirs are bound by the following 
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■'^PPei* and lowei* limits imposed "by the operating and 
engineering constraints. 

^i max max * (i=lf . - »p) » (k=l,...»N) (6.5) 

^i min -^i(k) max ' (i=P+l» • * • »P+li) , (k=l,...,N) 

( 6 . 6 ) 

^i min — ^ax ’ (i=l* • • . t p+h) , (k=2,...,N) 

(6.7) 

The formulation of the thermal subsystem and the electrical 
network is the same as in Sec. 2.5 given by e^ns.(2.15) to 
(2.26). Now the optimization problem is to find the optimal 
schedules for the pumped storage plants (pumping and genera- 
ting modes), conventional hydro plants and the thermal plants 
minimizing the fuel cost given by eqn.(2,15)» while satis- 
fying the various constraints (6.5) to (6.7)» (2.16), (2,17) 
and (2.20) to (2.26). 

6.3 DESCRIPTION OE THE ALGORITHM 

Step 1: To start with, a nominal hydro trajectory of 
storages of the pumped storage as well as the conventional 
hydro plants over the entire interval of optimization is 
assumed and is tested for feasibility as in Sec. 5.4. During 
the above procedure the discharge as well as pumping rates 
and their corresponding powers are available for each sub- 
interval for the nominal trajectory. 
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Step 2.: Considering the hydro generation or the pumping 
power required as injections or loads at the appropriate 
buses, optimal power flow solutions are obtained one for 
each subinterval to determine the corresponding optimal 
thermal generations in the system. The subinterval costs 
as well as total cost are calculated and stored. 

Step 3: The current nominal hydro trajectory is varied 
successively one hydro station after the other (including 
pumped storage plants) at all the instants k=2,...,lir 
using the ML7 as in Sec. 5*4 until a satisfactory conver- 
gence to the optimum solution is obtained. 

It is to be noted here that in the variational 
process, if at any instant k, the varied state value 
x^(k) > Xj^(k-l),the plant under consideration is in pumping 
mode, whereas if x^(k) < Xj^(k-l), it is in generating mode. 

6.4 FUMERICAl EESQITS 

The configuration of the system considered for 
illustration here is the same as in Figure 5.1 except that 
the hydro station at bus 1 is assumed to be a pumped 
storage plant. The interval of optimization (one day) is 
subdivided into 6 subintervals, each of 4 hours duration 
as in Sec. 5.5. The line and load data are the same as 
in Tables 5.1 and 5.2 respectively. Table 6.1 gives the 
hydro subsystem characteristics. The cost coefficients 
of the thermal stations are the same as in Table 5 ,4* 
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Table 6.2 gives the upper and lower limits on the reservoir 
storages at the various instants over the optimization 
interval. Table 6.5 stipulates the limits on the voltages, 
active and reactive powers at the various buses. Table 6.4 
gives the initial and optimal trajectories of the hydro 
reservoirs. Table 6.5 gives the optimal schedule for the 
pumped storage, conventional hydro and thermal plants 
obtained using the variable step size vector as in 
Chapters 4 and 5, with an initial step size of Ax = 200 
meter ^/sec.-4 hours for both the pumped storage and the 
conventional hydro plants. Table 6.6 gives the voltages, 
phase angles and reactive powers at all the buses over 
the entire interval of optimization (for the optimum 
condition). The pumping and generating efficiencies rip 
and T)„ are taken as 0.86 and 0.84 respectively. 



TABLB^ 6.1; CHABACTERISTICS OP THE HYDRO SHBHYSTEM 


Hydro system 

Pumped 
storage 
plant 
(i=l) ‘ 

Conventional 
hydro plant 

(i=2) 

Initial storage in meter 

sec. -4 hrs. 

2330 

2330 

Pinal storage in meter^/ 
sec. -4 hrs. 

2530 

2330 

Maximum allowable- discharge 
rate in met er^ /sec. 

500 

500 

Minimum allowable discharge 
ra.te in meter^/sec. 

-500 

100 

Initial head of the reser~ 



voir in meters 

14 

14 

Pinal head of the reservoir 

in meters 

14 

14 

Basic head 

01 

10 

10 

Inflow (assumed constant 
throughout) 

0 

300 


Rote: Negative sign in the above table shows the 
p-umping mole. 



TABLE 6.2: UPPER AND LOWER LIMITS ON THE RESERVOIR STORAGES 
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TABLE 6.3: LIMITS ON VOLTAGES AND GENERATIONS 
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TABLE 6.4; INITIAL AND OPTIMAL TRAJECTORIES OP THE RESERVOIR STORAGES 
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Note* The units are the same as defined earlier 



TABLE 6,5; OPTIMAL SCHEDULE FOR PUMPED STORAGE, HIDRO AHD THERI-dAL STATIONS 
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discharges Uj^(k) and u^Ck), 'which are given in meter ^/sec 

Initial cost in dollars = 11063.80 
Final cost in dollars = 11018,08 
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TABLE 6.6: VOLTAGES, PHASE ANGLES AND REACTIVE POWERS 
AT THE VARIOUS BUSES 


Snh- 

int erval 
k 

1 

2 

5 

4 

5 

6 

^1 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 


0.358 

-0.107 

0.583 

1.967 

0.339 

- 1.237 

^2 

o 

• 

H 

1.04 

1.04 

1.04 

1.04 

. 1.04 


-0.180 

-0.028 

-0.025 

-0.387 

1.399 

0.175 


1.0275 

1.0248 

1.0219 

1.0229 

1.0256 

1.0284 


-0.104 

-0.495 

-0.922 

-0.659 

1.182 

0.698 

u 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 


6.5 

5.991 

5.784 

7.378 

9.369 

8.64 

U 

1.04 

1.04 

1.04 

1.04 

1.04 

1.04 

s 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


-0.1995 

-0.1817 

-0.2097 

- 0.2703 

-0.1984 

-0.1219 

«2 

0.1386 

0,1465 

0.1636 

0.1815 

0.1187 

0.1378 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

«4 

-0.0747 

-0.0673 

-0.0616 

-0.0716 

-0.0793 

-0.0847 


0.2926 

0.2995 

0.3504 

0.4049 

0.3576 

0.2253 


Note: Q^, are the reactive powers at "buses 

respectively. All units are in p.u. except the phase 
angles (6’s), which are given in degrees. 
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6.5 Discussion on the series COnnECTED hydro reservoirs 

Case (i) Witho-ut Transport Delays: 

In all the problems considered hitherto, in this 
thesis, the hydro reservoirs are considered to he inde- 
pendent of each other. The extension of the solution 
methods discussed so far to the case of series connected 
reservoirs is straightforward as can he seen in this 
section. Consider a power system with h hydro plants with 
a series operation of some of the reservoirs. The dynamics 
of the hydro system may he represented as 


h 

Xjj^(k+1) = + D^(k) - u^(k) + I 

I 

( i— U- ,...,h), ( k=l , . . • ,n ) 


u. . u . (k) , 


( 6 . 8 ) 


where 


1 if reservoir i is directly upstream from 
reservoir D 

(6.9) 

0 otherwise 


The operating constraints remain the same as in the earlier 
cases as 


"^i min 

X . 

1 min 


< Uj_(k) < 

£ 3:^(k) < 


u. 

1 max 
^i max 


(i— 1, . . . ,h) , (k— 1, 
(i=l,.. .,h), (k=2. 


,H) (6..10) 

,IT) (6,11) 


The algorithms of Secs. 4. 7, 5*4 and 6.3 can he carried out 
with the following simple modification. In the variational 
process, let the state value corresponding to the 1-th 
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reservoir at the instant k he •under consideration. When 
x^(k) is incremental to Xj^(k) +Ax, the discharge rates 
in the (k--l)-thand k-th s-ubintervals are e'valuated using 
(6.8) for the i~th as well as the reservoirs situated 
downstream to the i-th reservoir. If any of the discharge 
rates or the storages of the reservoirs under consideration 
violate the limits on them, the variation under considera- 
tion is declared infeasible and the process 'is repeated 
as in the earlier cases. The rest of the steps in the 
above algo^rithms remain unaltered. 


Case (ii) With Transport Delays: 

Let T. be the transport delay (which is assumed to 
J 

be expressible in terms of integral multiple of the sub- 
interval duration), which elapses before the water released 
from the j-th reservoir reaches the i-th reservoir (down- 
stream of j-th reservoir). Then the dynamics of the hydro 
system may be represented as 


Xj_(k+1) = 2 :^(k) + Ij^(k) - u^(k) + j . u (k-v ) , 


0=1 3 

(i=l,...,h), (k=l,...,H) 


( 6 . 12 ) 


where is defined the same way as in (6.9). 

In applying the MIV in this case with transport 
delays, t:he following modification is to be made in the 
basic algorithms described in Secs. 4. 7» 5.4 and 6.3. 
Let the state value corresponding to the i—th reservoir 
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at the k~th instant he tinder consideration. ¥hen 
IS incremented to discharge rates in 

(k-l)-th and k-th subintervals along with discharge rates 
in the affected suhintervals of the reservoirs which are 
downstream to the i^th reservoir are determined using 
(6.11). The limits on the storages and discharge rates 
corresponding to these suhintervals are checked for the 
feasibility of the Variation under consideration. The 
costs of these suhintervals under consideration are 
determined (in the game way as in Secs. 4.7, 5.4 and 6.3) 
and com^pared with their costs before this variation is 
effected. Depending on the success or failure of this 
variation at_ the k-th instant, the process is advanced 
to the (k+l)-th instant. Rest of the steps remain the 
same as in the algorithms in Secs. 4.7, 5.4 and 6.3. 

6.6 CONCLUSIONS 

Numerical results are presented in Tables 6.1 to 
6.6^for a 5 bus system (whose configuration is the same 
as in Figure 5.1) consisting of one pumped storage, one 
conventional hydro and two thermal stations using the MLY 
and the Optimal Power piow method. The generality of the 
MLV in dealing with the pumped storage plants is established. 
Ill the engineering and operating limits on the hydro and 
pumped storage plants and the electrical network can be 
taken into account without any modification of the MLV 
and Optimal Power Plow algorithms discussed in Chapter 5. 
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A brief discussion on the consideration of the 
series connected reservoirs (without and with transport 
delays) in the optimal scheduling problem in hydro- thermal 
power systems is also included. 



CHAPTIR 7 


CORCIUSION 


7.1 INTRODUGTIOI 

This . chapt er is aimed at reviewing the significant 
results obtained during the course of this work and making 
a few suggestions for the future line of research in this 
area. Before reviewing the work done, a brief account of 
the objectives of this investigation is summarized. 

As is already mentioned in Chapter 1, the problem 
of optimal scheduling in hydro-thermal power system is 
more complex than in a purely thermal system, since the 
dynamics of the hydro reservoirs make the problem a 
variational problem, whereas a purely thermal scheduling 
problem is a static optimization problem. Several attempts 
have been made in the past to solve both the long range 
and short range optimal scheduling of hydro-thermal system 
problems using direct and indirect search methods. The 
application of direct search techniques like the Dynamic 
Programming [l5 ] » Incremental Dynamic Programming [lO] 
etc. have resulted in solution procedures requiring large 
computational time and .core storage. With larger inter- 
connected power systems, the procedures are hence 
hence inapplicable. Indirect search techniques using 



126 


the Pontriagin's maximum principle (both the continuous 
a,nd the discrete versions) result in a !Bwo Point Boundary 
Value Problem (PPBVP) to be solved iteratively. Phe lack 
of proper techniques for solving the above TPBVP satisfying 
the constraints on both the state and control variables 
make these methods applicable only to particular systems. 
Further the limits on the state variables present consi- 
derable difficulty in these methods. 

Hence the main objective of the present work has 
been to find a suitable solution technique which can be 
applied to the optimal scheduling problem in larger inter- 
connected hydro-thermal systems and which requires lesser 
computational time and core storage than the existing 
m.ethods, 

7.2 REVIEW OP THE WORK DONE 

The proposed solution method in this thesis 
consists of a decomposition of the combined systems into 
hydro and thermal subsystems and starting from a nominal 
hydro trajectory, the optimal trajectory is obtained by 
iteratively varying the initial nominal trajectory by a 
direct search technique known as MiV [23 ]» [24 ]. In each 
subinterval of optimization, a mathematical programming 
problem is to be solved for obtaining the optimal thermal 
generations. Augmented Penalty Punct ion Method [26 ] or 
Optimal Power Plow Method is employed to solve the above 
problem encountered, depending on the model 
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The first part of the thesis is devoted to investi- 
gate the solution using an approximate model discussed in 
Sec. 2.2. In this model the electrical network is repre- 
sented by the loss coefficients and the constraints on 
the voltage magnitudes, reactive powers and the line flows 
etc. have been ignored. To illustrate the method of 
solution a numerical example similar to the one considered 
in [lO] has been solved. In a on e-hydro- one- thermal problem 
considered here, the thermal generations are obtained by 
simply solving a quadratic equation. A detailed comparison 
of the solution algorithm based on MLV is made with that of 
IDP due to Bernholtz and G-raham [lO], The motivation for 
this comparison is that both these techniques are direct 
search techniques in the neighbourhood of a nominal 
trajectory. The superiority of the proposed algorithm over 
IIP procedure from the aspect of core storage requirements 
is established. Then the procedure is extended to a multi- 
hydro multi-thermal optimal scheduling problem, A numerical 
example of three hydro and four thermal stations has been 
solved to illustrate the solution procedure. In each 
subinterva,! the mathematical programming problem encountered » 
is solved by the Augmented Penalty Function Method to solve 
for the optimal thermal generations. With the simplicity 
and speed of this method, the overall solution procedure 
has become simple and fast. The problems solved in this 
part sho-w that the proposed method is quite promising 
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"based on the experience gained during the course of this 
work, it is recommended to start the MIY algorithm initially 
with as large a step size vector as is consistent with the 
problem and gradually reduce it each time as the iterative 
process converges. This procedure has produced significant 
reduction in computation times, compared to working with a 
constant (small) step size vector. The reduction is more 
significant in the larger dimensional problems (Secs* 4 . 8 , 
5.5) compared to the one-hydro-one-thermal case. 

To show the generality of the suggested solution 
techniques in the earlier parts, a general optimal schedule 
ing problem in a system consisting of pumped storage, 
conventional hydro and thermal plants is considered in 
Chapter 6 . The solution procedure does not need any 
modifications of the basic solution algorithms already 
employed in the earlier parts and the core storage and 
computational time requirements do not get altered 
significantly in this case. 

To summarize, the a.pplication of the MLY to 
correct the initial nominal trajectory iteratively 
in combination with either the Augmented Penalty Function 
Method or the Optimal Power Flow Method affords a simple, 
efficient and easily implemen table solution 'technique to 
the optimal scheduling problem in hydro- thermal power 
systems. 
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Based on the investigations made and the resxilts 
obtained, the following major conclusions can be drawn. 

(i) The algorithms proposed in Chapters 4 and 5» 

using the approximate and sophisticated models respectively^ 
provide simple and efficient procedures of solving the 
optimal scheduling problem for hydro-thermal power systems. 
In both the cases all the constraints are efficiently taken 
care of. Both the methods are applicable to problems of 
large interconnected systems, as the computational time 
and core storage requirements of these methods are signi-^ 
ficantly less compared to those of the existing methods 
[lO], [20] etc, 

(ii) The variable step size procedure suggested, whil© 
working with MIV, affords a significant reduction in the 
oomp-utational time required than that of working with 
constant (small) step size throughout. 

(iii) The inclusion of pumped storage plants in the 
proposed solution algorithms poses no difficulty. 

(iv) The extension of the solution algorithms 
sugg-ested need little modification to consider the 
se 3 ri. es connection of the hydro reservoirs with an 4 without 

tro-HSpci*"^ fime delays, 

7.5 IFUDUSE BINE OE RESEARCH 

As a sequel to the results obtained during the 

pres investigation, the following line of future 

resss3-3rch seems to be worth pursuing. 
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1. In all the problems considered, the load demand on 
the system and the inflows into the reservoirs have been 
taken to be deterministically known in advance. But in 
actual practice, these can only be characterized probabi- 
listically (particularly in the case of a long range 
scheduling problem) and hence the scheduling problem 
becomes a stochastic control problem. Suitable modifica- 
tions of the solution algorithms proposed in this thesis 
may be possible to consider the stochastic nature of the 
load demand and the inflow data. 

2. In the problems considered in Chapters 4 and 5, the 
objective function to be minimized and the equality con- 
straint were nonlinear and hence nonlinear programming 
techniques involving substantial computational effort 

had to be employed. Simple linear programming problems 
would result by linearizing the objective function and the 
equality constraint around the operating point, which can 
be solved more easily. Since working with the linearized 
equations will ensuire the satisfaction of the original 
constraints only in the neighbourhood of the operating 
point, reinitialization ha.s to be done once after every 
few iterations. An attempt by the author in this direction 
for a one-hydro-one- thermal problem seems to be quite 
promising and it appears that the linearization may lead 
to considerable reduction in computational times for 
higher dimensional problems. 
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